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Abstract 

In  geometrically  nonlinear  theories  of  inelasticity  of  solids,  the  deformation  gradient  is  typically  split  multiplicatively  into 
two  (or  more  terms),  none  of  which  need  be  integrable  to  a  motion  or  displacement  field.  Such  terms,  when  not 
integrable,  are  termed  anholonomic,  and  can  be  associated  with  intermediate  configuration(s)  of  a  deformed  material 
element.  In  this  work,  aspects  of  tensor  calculus  associated  with  anholonomic  deformation  are  analyzed  in  general 
curvilinear  coordinates.  Various  linear  connection  coefficients  for  intermediate  configurations  are  posited  or  derived; 
of  particular  interest  are  those  mapped  coefficients  corresponding  to  the  choice  of  identical  basis  vectors  in  multiple 
configurations.  It  is  shown  that  torsion  and  curvature  associated  with  such  mapped  coefficients  do  not  necessarily 
vanish,  even  though  torsion  and  curvature  tensors  of  the  original  connections  vanish  by  definition  in  reference  or 
current  configurations.  Intermediate  connection  coefficients  defined  in  this  way  exhibit  vanishing  covariant  derivatives 
of  corresponding  metric  tensors,  but  are  time  dependent  even  when  reference  (current)  configuration  connections  are 
fixed  in  time  at  a  given  material  (spatial)  location.  Formulae  are  derived  for  total  covariant  derivatives  of  two-  and  three- 
point  tensors  with  one  or  more  components  referred  to  the  intermediate  configuration.  It  is  shown  that  in  intermediate 
coordinates,  neither  the  divergence  of  the  curl  of  a  vector  field  nor  the  curl  of  the  gradient  of  a  scalar  field  need  always 
vanish.  The  balance  of  linear  momentum  for  a  hyperelastic-plastic  material  is  examined  in  the  context  of  curvilinear 
intermediate  coordinates. 


Keywords 

anholonomic  deformation,  curvature,  differential  geometry,  finite  strain,  multiplicative  decomposition,  torsion 


I.  Introduction 

The  notion  of  a  locally  relaxed  or  stress-free  intermediate  configuration  is  widely  used  in  geometrically  non¬ 
linear  (i.e.  finite  deformation)  models  of  solids.  An  intermediate  state  in  which  deformations  of  neighboring 
material  elements  may  be  incompatible  was  proposed  for  anelastic  materials  [1].  A  multiplicative  decomposi¬ 
tion  of  the  deformation  gradient  into  elastic  and  plastic  parts  was  developed  for  crystalline  solids  containing 
continuous  distributions  of  dislocations  [2,  3].  The  torsion  tensor  constructed  from  gradients  of  elastic  or  plastic 
deformations  can  be  associated  with  the  density  of  dislocations  [3-5].  Shortly  thereafter,  theories  of  continu¬ 
ous  bodies  with  inhomogeneities,  in  which  various  non-integrable  deformation  mappings  are  introduced,  were 
developed  [6,  7],  Subsequent  literature  on  nonlinear  kinematics  of  crystals  with  defect  distributions  is  immense; 
relevant  more  recent  works  include  [8-14]. 


Corresponding  author: 

JD  Clayton,  RDRL-WMP-B,  US  Army  Research  Laboratory,  Aberdeen  Proving  Ground,  MD  21005-5066,  USA. 
Email:  john.d. clayton  I  .civ@mail.mil 


Downloaded  from  mms.sagepub.com  at  US  ARMY  RESEARCH  LABORATORY  on  September  17,  2012 


Clayton 


703 


The  multiplicative  decomposition  of  the  deformation  gradient  for  crystals  deforming  via  dislocation- 
mediated  plasticity  has  been  extended  to  describe  kinematics  of  a  number  of  other  physical  phenomena.  For 
example,  such  decompositions  have  been  used  to  describe  thermal  expansion/contraction  [15],  growth  mechan¬ 
ics  in  biological  tissues  [16],  volume  changes  associated  with  point  defects  [10,  17]  and  voids  [18],  deformation 
twinning  [19,  20],  and  disclination  defects  [21,  22],  Theories  of  inhomogeneous  bodies  have  been  applied  to 
describe  solid  and  fluid  crystals  and  laminated  composite  materials  [7,  23]  in  addition  to  isotropic  solids  and 
anisotropic  solid  crystals.  Extensions  of  the  multiplicative  decomposition  to  three  (or  more)  terms  have  been 
used  to  quantity  elastic-plastic  phenomena  at  different  length  scales  [24],  Reviews  with  more  comprehensive 
lists  of  references  are  available  elsewhere  [14,  25].  It  should  be  noted  that  theories  for  nonlinear  material  behav¬ 
ior  in  the  context  of  defective  crystals  [8]  and  biological  systems  [26]  have  been  developed  that  do  not  require 
explicit  use  of  a  multiplicative  decomposition. 

In  the  present  work,  the  total  deformation  gradient  F  is  decomposed  as  F  =  FF,  where  F  and  F  are  generic 
two-point  mappings  from  the  reference  to  intermediate  configuration  and  intermediate  to  current  configuration, 
respectively.  In  coordinates,  =  d<pa/ dXA  is  the  tangent  map  from  the  reference  to  current  configuration, 
where  spatial  coordinates  xa  =  <pa(X,  i)  follow  the  motion  <p  that  may  depend  on  time  t  and  material  particle 
X.  Deformation  gradient  F  is  said  to  be  integrable  or  ‘holonomic’  since  the  differentiable  one-to-one  mapping 
<p(X,  t)  exists  between  referential  and  spatial  positions  of  material  particles.  On  the  other  hand,  F“A (X,  t)  and 
F~lo[a{x,  t)  need  not  be  integrable  functions  of  XA  and  xa,  respectively;  in  such  cases,  these  mappings  are  said 
to  be  ‘anholonomic’  [5,  14,  27,  28].  Correspondingly,  when  F  (F-1)  is  not  integrable,  coordinates  xa  that  are 
differentiable  one-to-one  functions  of  XA  (xa )  do  not  exist,  and  the  intermediate  configuration  is  said  to  be  an 
anholonomic  space. 

According  to  Schouten  [27],  the  study  of  differential  geometry  of  anholonomic  spaces  was  initiated  by 
Vranceanu  [29].  General  mathematical  identities  are  derived  or  listed  in  [27];  anholonomic  coordinate  transfor¬ 
mations  are  also  mentioned  by  Ericksen  [30,  p.  801]  in  the  context  of  tensor  fields  in  mechanics.  Applications 
of  anholonomic  geometry  to  plasticity  theory  are  described  by  Kondo  [5,  31], 

In  a  previous  work  co-authored  by  the  present  author  [28],  several  different  metric  tensors  on  the  intermediate 
configuration  were  discussed.  Of  particular  interest  in  that  work  was  a  metric  tensor  formed  from  the  scalar 
product  of  spatial  basis  vectors  convected  to  the  intermediate  configuration  using  the  elastic  distortion.  This 
metric,  which  corresponds  to  the  right  Cauchy-Green  deformation  tensor  formed  from  the  elastic  deformation 
gradient,  demonstrates  a  non-vanishing  curvature  tensor  associated  with  incompatibilty  or  anholonomicity  of 
the  intermediate  configuration,  and  proves  useful  in  formulating  scalar  energy  potentials  for  crystalline  solids 
with  continuous  distributions  of  dislocations  [14,  28].  In  that  work  [28],  the  same  notation  was  used  for  different 
sets  of  anholonomic  and  convected  basis  vectors;  the  present  work  corrects  this  ambiguous  notation. 

The  usual  choice  of  basis  vectors  in  the  intermediate  configuration  is  a  Cartesian  system  with  metric  ten¬ 
sor  components  equivalent  to  Kronecker’s  delta  symbols.  In  theories  of  crystal  elasto-plasticity,  coincident 
Cartesian  coordinates  are  almost  universally  used  for  every  configuration  of  a  material  body.  Cartesian  coordi¬ 
nates  prove  especially  convenient  for  representing  anisotropic  elastic  and  plastic  behavior  of  single  crystals,  for 
which  elastic  moduli  and  slip  system  geometry  are  most  easily  described  using  Cartesian  frames  of  reference 
[14,  25,  32,  33].  Teodosiu  [34]  suggested  an  intermediate  coordinate  system  wherein  an  orthonormal  director 
triad  is  attached  to  each  locally  relaxed  material  point,  but  triads  at  different  material  points  can  differ  by  a 
finite  rotation;  the  field  of  such  rotations  presumably  may  be  discontinuous  in  material  coordinates.  In  such  a 
representation,  the  metric  tensor  also  reduces  to  Kronecker’s  delta,  and  it  was  assumed  that  connection  coeffi¬ 
cients  associated  with  covariant  differentiation  in  the  (possibly  anholonomic)  intermediate  configuration  vanish 
identically  [34], 

Of  particular  interest  in  the  present  paper  are  general  curvilinear  coordinates  in  the  intermediate  configu¬ 
ration.  Curvilinear  coordinates  such  as  cylindrical  or  spherical  systems  are  useful  for  describing  bodies  whose 
shapes  are  naturally  extrinsically  or  intrinsically  curved,  e.g.  cylinders,  spheres,  and  shells  of  various  kinds. 
Such  shapes  arise  frequently  in  biological  systems.  Furthermore,  certain  crystalline  structures  may  be  amenable 
to  description  with  curvilinear  coordinates;  for  example  screw  dislocations  perpendicular  to  the  basal  plane 
in  crystals  with  hexagonal  symmetry  (with  second-order  elastic  constants  having  transverse  isotropy)  may  be 
described  using  cylindrical  coordinates  [35],  In  such  situations,  it  is  natural  for  one  to  use  the  same  system  of 
coordinates  in  all  configurations  if  the  shape  of  the  body  remains  roughly  the  same  during  deformation.  For 
example,  if  the  body  remains  cylindrical  throughout  the  deformation  process,  the  natural  choice  of  coordinates 
would  be  cylindrical  coordinates  in  reference,  intermediate,  and  current  configurations.  Simo  [36]  developed  a 
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finite  elastic-plastic  theory  wherein  the  same  metric  tensor  and  presumably  the  same  coordinate  system,  which 
may  generally  be  curvilinear,  are  used  for  both  reference  and  intermediate  configurations.  The  present  author 
[14,  p.  91]  suggested  that  two  pragmatic  choices  for  intermediate  coordinate  systems  are  either  this  prescrip¬ 
tion  [36]  or  the  prescription  of  the  same,  possibly  curvilinear,  coordinate  system  in  intermediate  and  spatial 
configurations. 

The  present  paper  supplements,  refines,  and  substantially  extends  a  previous  treatment  of  anholonomic 
geometry  in  the  context  of  solid  mechanics  (in  particular,  brief  sections  2.8  and  3.2.3  of  [14]).  Choices  of 
coordinate  systems  (i.e.  basis  vectors)  in  the  potentially  anholonomic  intermediate  configuration  are  critically 
examined.  Metric  tensors,  connection  coefficients,  torsion,  and  curvature  are  derived  for  each  choice  of  basis. 

As  shown  for  what  appears  to  be  the  first  time  in  this  paper,  the  choice  of  such  mapped  intermediate 
coordinate  systems  from  curvilinear  referential  or  spatial  coordinate  systems  [14,  36]  leads  to  intermediate 
basis  vectors  whose  derived  connection  coefficients  may  be  non-symmetric  in  covariant  indices  and  may  have 
non-vanishing  torsion  and  curvature.  It  is  also  shown  that  torsion  and  curvature  may  be  non-zero  even  if  the 
intermediate  configuration  is  holonomic.  Upon  development  of  logical  definitions  for  partial  and  covariant  dif¬ 
ferentiation  with  respect  to  possibly  anholonomic  intermediate  coordinates,  it  is  shown  that  the  divergence  of 
the  curl  of  a  vector  field  and  the  curl  of  the  gradient  of  a  scalar  field  need  not  necessarily  vanish  in  the  interme¬ 
diate  configuration.  Piola’s  identity  for  the  Jacobian  determinant  [14,  37,  38]  also  does  not  generally  apply  in 
the  intermediate  configuration.  Formulae  are  derived  for  total  covariant  derivatives  [30,  39]  of  two-  and  three- 
point  tensors  with  one  or  more  components  referred  to  the  intermediate  configuration.  It  is  shown  by  example 
that  such  formulae  are  needed  when  writing  the  balance  of  linear  momentum  for  nonlinear  hyperelastic-plastic 
solids  in  general  curvilinear  coordinates. 

Much  of  this  paper  is  tutorial  in  nature  (e.g.  the  content  of  Section  2  can  be  found  in  books  on  nonlin¬ 
ear  continuum  mechanics  and  elasticity  in  curvilinear  coordinates);  however,  in  addition  to  serving  as  a  useful 
reference,  such  content  is  needed  to  develop  and  contrast  results  in  subsequent  sections.  As  mentioned  in  the 
preceding  paragraph,  the  present  work  contains  apparently  new  results  pertaining  to  particular  choices  of  inter¬ 
mediate  coordinate  systems  that  are  not  explicitly  evident  in  other  relevant  works  incorporating  direct  (i.e. 
coordinate-free)  notation  [6,  40], 

This  article  is  organized  as  follows.  Definitions  and  notation  for  geometry  and  kinematics  of  holonomic 
(i.e.  integrable)  deformation  gradients  are  given  in  Section  2.  Kinematics  and  geometry  of  anholonomic  defor¬ 
mation  are  described  in  Section  3,  including  integrability  conditions,  choices  of  coordinate  systems,  rules  for 
partial  and  covariant  differentiation,  and  important  derived  identities.  The  local  balance  of  linear  momentum  in 
geometrically  nonlinear  continuum  mechanics  is  examined  in  Section  4  in  the  context  of  possibly  curvilinear 
intermediate  coordinate  systems.  Conclusions  follow  in  Section  5,  while  an  Appendix  contains  explicit  forms 
of  intermediate  connection  coefficients  and  curvature  for  cylindrical  coordinates.  Notation  of  nonlinear  contin¬ 
uum  mechanics  [14,  38,  39]  is  used.  Einstein’s  summation  convention  applies  over  repeated  indices.  Uppercase 
Roman  font  is  used  for  indices  corresponding  to  referential  (i.e.  material)  coordinates,  lowercase  Roman  font 
for  current  (i.e.,  spatial)  coordinates,  and  lowercase  Greek  font  for  intermediate  (and  possibly  anholonomic) 
coordinates. 

2.  Holonomic  deformation 

2.1.  Configurations,  coordinates,  and  metrics 

A  material  point  or  particle  in  reference  configuration  B0  is  labeled  X.  The  corresponding  point  in  the  current 
or  spatial  configuration  B  is  labeled  x.  Time  is  denoted  by  t.  The  motion  and  its  inverse,  respectively,  are 

xa  =  <pa(X,  t)  =  xa(X,  0,  XA  =  O  'U,  t )  =  XA(x,  t).  (1) 

Unless  noted  otherwise,  spatial  coordinates  xa  and  reference  coordinates  XA  are  assumed  sufficiently  differen¬ 
tiable  functions  of  their  arguments.  Partial  differentiation  is  written  alternatively  as  follows: 


3(0  a  n  r, 

—  °aC)  — 

3(0  «  M  , , 

dx“  ~  afl°  ~ 

(2) 

The  following  identities  are  used  frequently: 

3^[3s(-)]  =  35[3^(-)], 

3fl[3fe(0]  =  3*[3a(-)]. 

(3) 
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Let  X  e  B0  and  x  e  B  denote  position  vectors  in  Euclidean  space.  Referential  and  spatial  basis  vectors  are  then 
written,  respectively,  as 

G A(X)  =  dAX,  g  a(x)  =  3flx.  (4) 

Reciprocal  basis  vectors  are  GA(X)  and  ga(x).  Scalar  products  of  basis  vectors  and  their  reciprocals  are 

(Ga  ,Gb)  —  <gfl,g  b)=8ab.  (5) 

Kronecker  delta  symbols  are  <5#  and  8ah.  Indices  in  parentheses  are  symmetric,  indices  in  square  brackets  are 
anti-symmetric,  and  indices  between  vertical  bars  are  excluded  from  (anti-)symmetry  operations,  e.g. 


A(ab)  —  ~(^ab  +ABa),  A[ab\  —  - (Aab  —^ba),  A^\c\b]  —  ~(AAcb  ~  ABca)- 


From  (3), 


9[^[3^](-)]  =  3[a[3*](-)]  =  0. 


(6) 

(7) 


Using  (5),  the  scalar  product  of  a  vector  V  =  VAGA  and  a  covector  a  =  aB GB  is 

<V,a>  =  VAaB(GA,  GB)  =  VAaB8B  =  VAaA.  (8) 

The  tensor  product  or  outer  product  ®  obeys 

(G4  ®  G5)Gc  =  G^IG5,  Gc>.  (9) 


Symmetric  metric  tensors  G(X)  and  g(x)  are  introduced  for  respective  configurations  B0  and  B: 

G(X)  =  GabGa®Gb  =  (Ga-Gb)Ga®Gb,  GAB  =  G(AB);  (10) 

g(x)  =  gabg“  ®  gb  =  (ga  '  gb)gC'  ®  g\  gab  =  g(ab),  (1  1) 

where  the  dot  product  of  vectors  V  and  W  is 

V-W=VA  Wb(Ga  ■  GB)  =VaWa.  (12) 

As  indicated,  metric  tensors  can  be  used  to  lower  contravariant  indices: 

VA  =  VbGab,  Ga  =  GabGb.  (13) 

Inverses  of  metric  tensors  are 

G-1  =  GabGa  ®Gb  =  (G"4  •  Gb)Ga  ®  Gb,  GAB  =  G(AB)\  (14) 

g”‘  =  g“bga  ®gb  =  (gfl  •  gb)ga  ®  gb,  gab  =  g{ab)\  (15) 

where  the  dot  product  of  covectors  a  and  fi  is 

ct-P  =  aAMGA-GB)  =  aA/3A.  (16) 

As  indicated,  inverse  metric  tensors  can  be  used  to  raise  covariant  indices: 

VA  =  VbGab,  Ga  =  GabGb.  (17) 

Furthermore,  from  the  definition  of  the  inverse  operation, 

GACGCB  =  8j,  gacgcb  =  8ab.  (18) 
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Determinants  of  metric  tensors  and  their  inverses  are  written  as 

G  =  det  G  =  det  (GAB),  G~l  =  det  G1  =  det  (GAB)  =  1/G;  (19) 

g  =  det  g  =  det  ( gab ),  g~ 1  =  det  g“ 1  =  det  (gah)  =  1  /g.  (20) 

Permutation  tensors  are  defined  as 

eABC  _  gABC ,  =  J~G  eABC',  ^  ^ ^ abc  —  y/g  C abc ■  (21) 

VG  VS 

Permutation  symbols  are  eABC,  eABC,  eabc,  and  eah(:;  definitions  and  identities  involving  third-order  permutation 
symbols  implicitly  assume  a  three  dimensional  ambient  space.  In  this  paper,  metric  tensors  are  assumed  positive 
definite  over  any  volume  [39];  however,  in  certain  coordinate  systems,  determinants  of  metric  tensors  or  their 
inverses  may  be  zero  or  undefined  at  certain  points,  lines,  or  surfaces;  for  example,  along  null  radial  coordi¬ 
nates  in  spherical  or  cylindrical  systems.  Shifter  tensors,  examples  of  two-point  tensors,  can  be  introduced  in 


Euclidean  space  [30,  39]: 

gaAV,X)=(ga,GA),  gAa(x,X)=  {ga,GA);  (22) 

gaA(x,X )  =  ga  ■  Ga  =  gia,  gaA(x,X)  =  ga  ■  GA  =  gAa;  (23) 

gsgb  =  K’  SbglB  =  ’  (24) 

gaA  =  gabgA  =  G/iBg^  =  gabGABgBb',  (25) 

det  (gf)  =  1  /  det  (gaA)  =  V det  (gab)/  det  (GAB)  =  VgjG.  (26) 

Vector  V  and  covector  a  can  then  be  written 

V  =  VaGa  =  (VaVa)(gbAgb)  =  Va8bagh  =  Vaga;  (27) 

a  =  aAGA  =  (aagA)(gt  gb)  =  oia8abgb  =  aag“.  (28) 

The  following  rules  for  shifting  of  basis  vectors  are  implied: 

g  a  =  gfGA,  ga  —  gAGA;  GA=gaAga,  Ga  =  gAgf‘.  (29) 


2.2.  Linear  connections 

Discussed  in  what  follows  are  connection  coefficients  and  covariant  differentiation,  torsion  and  curvature 
tensors,  Riemannian  geometry,  and  features  of  Euclidean  space. 


2.2.1.  Covariant  derivatives.  Associated  with  a  linear  connection  in  configuration  B0  is  the  covariant  derivative 
operator  V.  The  covariant  derivative  assigns  to  two  vector  fields  V,  W  a  third  vector  field  VvW,  called  the 
covariant  derivative  of  W  along  V.  In  coordinates, 

Vv  W  =  (VBdBWA  +  T  Vc  Wc  Vb)Ga  .  (3  0) 


In  ^-dimensional  space,  w3  entries  of  the  object  V'V  are  called  connection  coefficients.  The  covariant  derivative 
is  applied  to  components  of  vectors  and  tensors  of  higher  order  as  [27] 


VnAg;;m  —  9iv^G...M  +  r nrAg...m  + 


A...F 


.Jt  aR...F 


,  aa-r  _  a 
^  1  NR^G-.M  1  "r- 


,r4a...f 

ngar...m 


p..fl  aA...F 

1  nmag...r- 


(31) 


Indices  of  covariant  differentiation  on  the  left-hand  side  correspond  to  those  of  partial  differentiation  and  the 
first  covariant  index  of  the  connection  coefficients  on  the  right-hand  side.  The  covariant  derivative  of  a  scalar 
field  A (X )  is  equivalent  to  its  partial  derivative,  i.e.  VNA  =  dNA .  Definitions  can  be  applied  analogously  to 
configuration  B  with  coordinates  xa. 
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2.2.2.  Torsion  and  curvature.  The  torsion  tensor  of  an  arbitrary  connection  with  coefficients  T^i  is  defined  in 
holonomic  coordinates  as 

T  =  Tjg Gb  0  Gc  0  Ga  =  r^c]GB  ®GC®  G^.  (32) 

Components  of  the  torsion  tensor  are  the  anti-symmetric  covariant  components  of  its  corresponding  connection: 
=  r^q.  Components  of  the  Riemann-Christoffel  curvature  tensor  of  this  arbitrary  connection  are 


nmmmA  <V  yn . mA  i  rj  T~^  mmA  ri  m.E 

kBCD  —  za[Bl  C]D  At  [B\E\l  C]D- 


(33) 


Both  the  torsion  and  curvature  transform  under  a  conventional  change  of  basis  as  true  tensors.  Different  defini¬ 
tions  are  used  for  torsion  and  curvature;  those  listed  here  follow  [14,  27].  Skew  second  covariant  derivatives  of 
a  contravariant  vector  V  and  a  covector  a  can  be  expressed  as  [27] 


VfflVq VA  =  V[B(dC]VA  +  F$dVd)  =  l-R^DVD  -  T$VdVa,  (34) 

=  X[B(dc]CiD  —  Tq^a^)  =  —  -RBcD&a  ~  T'bc^ 'a<xd-  (35) 


2.2J.  Arbitrary  holonomic  connections.  Coefficients  of  an  arbitrary  connection  in  holonomic  coordinates  can  be 
written  [27] 


-ImmA 

BC 


=  {. 


.mA  1 


1  I  rJ~'**A  rriuA,  ,  m 

BC*  .BC 


l 

2( 


+  ^  Wc  +  -  M'Bc). 


(36) 


Here,  Gab  and  its  inverse  GAB  are  symmetric,  three  times  differentiable,  invertible,  but  otherwise  arbitrary 
second-order  tensors.  Christoffel  symbols  of  the  tensor  Gab  are 


{ )}A: }  =  -GAD(dBGcD  +  dcGBD  ~  <)[)Gbc)  (symbols  of  the  second  kind);  (37) 

Ibcj]  =  Gad{'bc}  (symbols  of  the  first  kind).  (38) 

The  third-order  object 

MSc  =  GadMbcd  =  -GadVbGcd  =  GCDVBGiD,  (39) 


where  the  final  equality  follows  from  V b(GqdGad)  =  VgSjA  =  0.  The  covariant  derivative  of  Gab  follows  as 


V ’aGBc  =  3aGBc  —  r aBGdc  —  T acGbd  =  —Mabc  =  —Ma(Bc). 


(40) 


When  VaGBc  =  0  (or  when  MABc  =  0),  the  connection  is  said  to  be  metric  with  respect  to  Gab,  he.  a  metric 
connection.  For  a  metric  connection,  covariant  differentiation  via  V  and  lowering  (raising)  indices  by  Gab  (GAB) 
commute. 


2.2.4.  Riemannian  geometry.  In  Riemannian  geometry,  by  definition,  the  torsion  vanishes  and  the  connection  is 
metric: 

hflc  =  Iflc)  =  ]2Gad(8bGcd  +  dcGBo  ~  BoGbc),  T$  =  0,  MABC  =  0.  (41) 

In  Riemannian  geometry,  the  number  of  independent  components  of  the  Riemann-Christoffel  curvature  tensor 
is  n2(n2  —  1)/ 12;  for  example,  one  component  for  two-dimensional  space  and  six  for  three-dimensional  space. 

2.2.5.  Euclidean  space.  Let  Gab(X)  =  Ga  ■  Gb  be  the  metric  tensor  of  the  space.  The  Levi-Civita  connection 
coefficients  of  Gab,  written  as  p  '^Ac,  are  the  associated  metric  and  torsion-free  connection  coefficients  of  (41): 

p  «c  =  \gAD(4bGCd  +  3 cGbd  -  3 DGBC)  =  f  fB.  (42) 
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The  superposed  G  is  a  descriptive  label  rather  than  a  free  index  and  is  exempt  from  the  summation  convention.  In 
Euclidean  space,  the  Riemann-Christoffel  curvature  tensor  of  the  Levi-Civita  connection  vanishes  identically: 

G  ..~4  9  /  a  G  .Jt  I  G  .Jt  G  ..E  \  n 

R  BCD  ~  Z  l  °[5p  C]D  p  [B\E\  f  C]D )  ~  U- 

In  ^-dimensional  Euclidean  space,  a  transformation  to  a  ^-dimensional  Cartesian  coordinate  system  is  admitted 
a  each  point  X  such  that  GAB(X)  — >  SAB,  where  SAB  are  covariant  Kronecker  delta  symbols.  When  the  cur¬ 
vature  tensor  of  the  connection  vanishes,  the  space  is  said  to  be  intrinsically  flat;  otherwise,  the  space  is  said 
to  be  intrinsically  curved.  In  two  dimensions,  a  cylindrical  shell  is  intrinsically  flat,  while  a  spherical  shell  is 
intrinsically  curved. 

Henceforward,  reference  configuration  B0  is  treated  as  a  three-dimensional  Euclidean  space.  More  precisely, 
in  this  manuscript,  when  a  configuration  is  identified  with  ^-dimensional  Euclidean  space,  a  simply  connected 
deformable  body  of  finite  size  in  that  configuration  is  assumed  to  occupy  an  open  region  of  infinitely  extended  n- 
dimensional  Euclidean  vector  space.  In  the  interest  of  brevity,  such  a  configuration  is  simply  labeled  a  Euclidean 
space. 

The  covariant  derivative  associated  with  the  Levi-Civita  connection  on  B0  is  written  as 

y/)  =  (•);,■  (44) 

Covariant  derivatives  of  basis  vectors  and  their  reciprocals  vanish,  leading  to 

Ga;b  =  dBGA  -  G  bcaGc  =  0  =>  dBGA  =  G  fA Gc,  (45) 

Ga.r  =  dBGA  +  cf  fcGc  =  0^dBGA  =  -^  fcGc.  (46) 

It  follows  that  Christoffel  symbols  can  be  computed  as 

f  'bc  =  r  BC^D  =  r  bc(Gd,  Ga)  =  ( dBGc ,  G^>.  (47) 

The  partial  derivative  of  the  metric  tensor  is 

dcGAB  =  dc(GA  ■GB)  =  2GTf(A GB)D.  (48) 

Because  the  Levi-Civita  connection  is  a  metric  connection, 

GAB-f  =  dcGAB  —  p  'caGdb  —  p  "cbGad  =  0.  (49) 

From  the  symmetry  of  the  Levi-Civita  connection  [or  (4)  and  (7)],  skew  partial  derivatives  of  basis  vectors 
vanish: 

3[itG^]  =  p  -{bca]Gc  =  d[BdA]X  =  0.  (50) 


The  following  identity  will  be  used  for  the  derivative  of  the  determinant  of  a  non-singular  but  otherwise  arbitrary 
second-order  tensor  A  [30]: 


3  det  A 

3  K 


=  A  det  A. 


(51) 


Applying  this  identity  to  the  determinant  of  the  metric  tensor, 


dA  (In  VG)  —  y  ba- 


(52) 


Now  consider  the  spatial  configuration  5,  which  is  also  henceforward  identified  with  Euclidean  space,  with 
metric  tensor  components  gab(x )  =  go  ■  g b-  The  spatial  Levi-Civita  connection  is 


g  ..a 

Y  be 


id^bgcd  T  3 egbd  (fig he)  — 


g  ..a 

P  cb  ■ 


(53) 
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Superposed  g  is  a  descriptive  label  rather  than  a  free  index  and  is  exempt  from  the  summation  convention.  In 
Euclidean  space,  the  Riemann-Christoffel  curvature  tensor  vanishes: 

g  ...a  o  (  a  g  ..a  ,  g  ..a  g  ..e  \  a  /c,r\ 

Jibed  —  Z  l  0P>r  c]d  m  Y  [b\e\Y  c]d )  ~  U- 

The  covariant  derivative  associated  with  the  Levi-Civita  connection  on  B  is  written  as 

W')  =  (■);-  (55) 


Analogously  to  identities  for  the  reference  configuration, 

ga;  b  =  0  =»  dhga  =  kr  Bcagc,  (56) 

8%  =  0  =>  hga  =  ~y  bcSc,  (57) 

Ybc  =  (dbgc,ga),  (58) 

=  r  E bdfic  =  =  o,  (59) 

degab  =  28rc(agb)d,  (60) 

9fl(lnv@  =  f^,  (61) 

gab\c  =  0.  (62) 


2.3.  Differentiable  operators 

Differential  operators,  specifically  the  gradient,  divergence,  curl,  and  Laplacian  are  defined  for  holonomic 
coordinates  in  Euclidean  space.  Partial  and  total  covariant  derivatives  of  two-point  tensors  are  introduced. 

2.3. 1 .  Gradient,  divergence,  curl,  and  Laplacian.  The  description  that  follows  in  this  subsection  is  framed  in  Euclidean 
reference  configuration  B0  spanned  by  holonomic  coordinates  XA.  Analogous  definitions  and  identities  apply 
for  a  description  in  the  spatial  configuration.  The  gradient  of  a  scalar  function  f(X)  is  equivalent  to  its  partial 
derivative: 

y/  =  y  JGA  =/,aGa  =  d^G4.  (63) 

The  gradient  of  a  vector  field  V(X)  =  V4  G/(  is 

c,  V  =  dBV  <g>  Gs  =  fdB  VA  +  °  <g>  Gb  =  VfB Ga  <S>  GB  =  y  4  VbGa  <g>  GB.  (64) 

The  trace  of  a  second-order  tensor  A  is 


trA  =  Aaa  =  AabGab  =  AabGab  . 


(65) 


The  divergence  of  a  vector  field  V(X )  =  VAGA  is 

(y,  V)  =  tr  (y  V)  =  (dAV,GA)  =  V^  =  dAVA  +  ^-/bVb  =  -j=dA(VGVA),  (66) 

where  the  final  equality  follows  from  (52).  The  vector  cross  product  x  obeys,  for  two  vectors  V  and  W  and  two 
covectors  a  and  /3, 

V  x  W  =  eABCVBWcGA,  a  x  ft  =  eABCaBpcGA.  (67) 

The  curl  of  a  covariant  vector  field  a(X )  obeys 


x«-GJx  dA(aBGB)  —  GA  x  GBaB-A  —  eABLac;BGA  —  eABL dBacGA 


ABC  ~ 


r ABC ; 


(68) 
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where  the  final  equality  follows  from  symmetry  of  the  Levi-Civita  connection.  The  Laplacian  of  a  scalar  field 
f(X)  is 

2  ^,AR  ^  ...  1  ,,  ,  rXt^AR^ 

(69) 


v  /  =  GABf-AB  =  {GabU\,b  =  —j=  dB(VGGAB djf). 


The  divergence  of  the  curl  of  a  (co)vector  field  vanishes  identically: 

/G  G  v  \  _  ABC  _  C\AB\  _  a 

\V’Vxa'  —  €  <*c;ba  —  g  aC;(AB)  —  v, 
as  does  the  curl  of  the  gradient  of  a  scalar  field: 


G 

V  ~  V 


x  V-/  —  CBCf-^cB^A  —  f'4LtfLJ/;(5C)G^  —  0. 


A[BC\{ 


(70) 


(71) 


2.3.2.  Partial  and  total  covariant  derivatives.  Consider  a  two-point  tensor  (i.e.  double  tensor)  A(X,x)  of  order  two: 

A  =  A^ga®GA.  (72) 

Components  of  the  total  covariant  derivative  of  A^  are  defined  as  [14,  30,  39] 

dxb 


(dAa 
~dX^ 


G  ..C  aci  |  , 

Y  baa.c  I  ' 


dxb 


I  S  ..a  ac  _ 

x+YbcAA  j  dxB- 


(73) 


Quantities  'dBxh  will  be  identified  with  components  of  the  deformation  gradient  in  Section  2.4.1.  Partial  covariant 
derivatives  of  A"a  are  defined  as  usual  covariant  derivatives  with  respect  to  indices  in  one  configuration,  with 
those  of  the  other  configuration  held  fixed: 


4a  _  ^ 
aa:b  ~ 


dXB 


G  ..C  a  a 
Y  BAA.C’ 


A« 

AJl\b  — 


dxb 


x 


+ 


g  ..a  ac 

r  bcAA 


From  (1),  xa  =  xa(X,  t).  Writing  A°4[X,x(X,  t)}  =  A^(X,  t),  the  partial  derivative  of /)",  at  fixed  t  is 


dXB  dXB 


dXB 


+ 


9^ 


dxh 


dBC 


x 


With  this  more  conventional  notation,  (73)  becomes 

dA“ 


« 4 


Aa  — 

■ A'B  dXB 


-  Aa  i  s  .. a  ac  a  0 

r  BAA.C  "t  p  bcAAaBx  ■ 

The  total  covariant  derivative  can  also  be  obtained  directly  by  inspection  of  the  (material)  gradient  of  A: 


yA  =  dBA  ®GB  =  dB(A^ga  ®  Ga)  ®Gb=  A^:b ga  ®ga®  gb. 
The  total  covariant  derivative  can  be  extended  to  two-point  tensors  of  arbitrary  order  as  [30] 

Gl'.'D'.x  = 


(74) 


(75) 


(76) 


(77) 


(78) 


where  partial  covariant  derivatives  are  taken  with  respect  to  indices  referred  a  single  configuration,  as  in  (74). 
From  the  definition  of  the  total  covariant  derivative,  noting  that  dKxk  djXK  =  Sf,  total  covariant  derivatives  map 
between  configurations  like  partial  derivatives: 


GaG):k\Xk  =  (AaA;fE);K  dkXK  +  (A^hdKX1  dkXK  =  (AaA;;eE):k. 


(79) 


2.4.  Kinematics  ofintegrable  deformation 

The  deformation  gradient  and  its  inverse  are  introduced.  The  Jacobian  determinant  associated  with  volume 
changes  is  defined,  and  Piola’s  identities  are  listed. 
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2.4.1.  Deformation  gradient.  The  deformation  gradient  F  is  the  two-point  tensor  with  components,  from  (1), 

F  =  F\ ga  <g>  Ga,  (80) 


„  dtp  (X,  t)  dxa(X,t ) 

FaA(x,  t)  =  y  =  — ~y  =  xaA  =  dAxa 
>  dxA  dxA  ,A 

Similarly,  the  inverse  deformation  gradient  and  its  components  are 

F  =  F~u  Ga  <g>  ga, 


F~lA.a(x,t)  = 


d&A(x,t)  dXA(x,  t) 


dxa 


dxa 


=  X/‘  =  daX  . 


(81) 

(82) 

(83) 


From  differentiability  and  invertibility  properties  of  the  mappings  in  (1),  detF  A  0  and  detF  1  A  0.  By 
definition, 

(84) 


i?a  r'—^A  _  oa  T7~^A  j?a  _  ?A 

r  .b  ~  f  .cr  ,B  ~  °B- 


Partial  differentiation  (holding  time  t  fixed)  proceeds  as 


dA(0  =  U-)dAxa  =  9flQ  F\,  da(-)  =  dA(-)daXA  =  dA(-)F~1Aa. 


(85) 


Consider  a  differential  line  element  dX  in  the  reference  configuration.  Such  an  element  is  mapped  to  its 
representation  in  the  current  configuration  dx  via  the  Taylor  series  [14,  41] 


dx\X )  =  {FU)\xdXA  +  yFy)\xdXA  dXB  +  ^{F:a.bc) 


dXA  dXB  dXc  +  • 


where  components  of  the  total  covariant  derivative  of  F  are  defined  as  in  (76): 


dFa 

Z7<z  _  AG  m.C  E1  a  |  g  .. a  j?c  j?b 

r  A'.B  ~  d  TRAr.r^~ 


dXB  rM'-cT  T  M  —  9 B(dAx‘  )  p  ba^cx  +  p  be  dAx  dgx  —  x:AB. 


Third-order  position  gradient  follows  likewise  as 

F^4‘bc  =  (F“a-b):c  =  (x- Ab):c  —  X-ABC  =  dc[dB(dAx  )]  + 


(86) 


(87) 


(88) 


From  the  identity  9^[9g(-)]  =  9«[9.i(-)]  of  (3)  and  symmetry  of  the  (torsion- free)  Levi-Civita  connection 
coefficients  in  both  reference  and  current  configurations,  it  follows  that 


r?a  _  r?a  _  rpa  r?a  _  rpa 

r  A'.B  ~  r  .B'.A  ~  r  .(A'.B)’  r  A'.BC  ~  r .( A'. 


To  first  order  in  dX,  (86)  is 


dxa(X)  =  (dAxa)  dXA  =  (FaJ 


x 


x 


.(a:bo- 


dXA,  dx  =  FdX. 


(89) 


(90) 


Relationship  (90)  is  the  usual  assumption  in  classical  continuum  field  theories  [37]. 

2.4.2.  Jacobian  determinant.  Jacobian  determinant  J  relates  differential  volume  elements  in  reference  and  current 
configurations: 

dv  =  JdV.  (91) 

Differential  volume  elements  in  the  reference  (dV)  and  current  configuration  (dv)  are  written  symbolically  as 

dV  =  yfGdX 1  dX2  dX 3  c  B0,  dv  =  gdx 1  dx2  dx 3  C  B.  (92) 
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The  Jacobian  determinant  J\¥(X,  t),g(x),  G(X)]  is,  from  (21)  [14,  37-39] 


T  —  -<ZABC<Z 

•J  —  ,  £  ^abc 


=  \VslG  eABCeahcF:AF!;BFfc  =  ^Jg  detF  =  / 


Aet(dAxa).  (93) 


From  similar  arguments,  inverse  Jacobian  determinant  J  1  [F  1  (x,  t),  G(X),  g(x)]  =  1  /./  is 
J~l  =  \eabceABC  F-lAMF-x\F~lcmC  =  ^det  (GAB)/  det  (gah)  det(3flX^). 


(94) 


When  motion  is  restricted  to  rigid  translation  (or  no  motion  at  all),  then  F"A  =  g"  [the  shifter  of  (22)],  F  = 
gAga  ®  GA,  and  J  =  y/g/G  dctfg")  =  1  follows  from  (26).  Since  volume  remains  positive,  J  >  0  and  det  F  >  0. 
From  (51), 

dJ  _ lA  dJ~x 


9^ 


=  JF 


dF 


—  T~lFa 
-IA  r-4' 


(95) 


One  of  Piola’s  identities  is  derived  by  taking  the  divergence  of  the  second  expression  of  (95)  [14]: 

=  da(J~lFaJ  +  J-lgr  -?bF\  -  J-xgt£  =  J~XF~lBma  [dB(dAxa)  -  dA(dBxa)-\  =  0. 

The  same  steps  performed  on  (95)  with  reference  and  spatial  coordinates  interchanged  produce 
(JF~lAMyA  =  JFbA  [dbF~lAa  -  daF~XAb\  =  JFb4  [dh(daXA)  -  3fl(3^)]  =  0. 

Let  the  vector  field  A(X)  =  AaGa  be  the  Piola  transform  of  a(x)  =  aa ga: 

Aa  =  JF~lAaaa. 

Taking  the  divergence  of  (98)  and  applying  the  product  rule  for  covariant  differentiation  along  with  (97)  gives 
4-a  =  A?a  =  (JF-lAJ:Aa“  +JF~XAaa“A  =  JF~IAXa  =  Ja°a  =  Ja“  (99) 


(96) 


(97) 


(98) 


3.  Anholonomic  deformation 

3.1.  Anholonomic  spaces  and  geometric  interpretation 

The  deformation  gradient  is  split  multiplicatively  into  two  terms,  neither  of  which  is  necessarily  integrable. 
Geometric  consequences  of  such  a  construction  are  considered,  including  rules  for  differentiation,  coordi¬ 
nate  systems,  metric  tensors,  and  connection  coefficients  convected  from  reference  or  current  configurations 
to  anholonomic  space. 


3.1.1.  Two-term  decomposition  of  deformation  gradient.  Consider  a  multiplicative  split  of  the  (total)  deformation 
gradient  F  into  two  terms: 

F  =  FF,  (100) 

or  in  indicial  notation,  letting  the  Greek  index  a  =  !,...«,  where  n  is  the  dimension  of  Euclidean  space, 


dcpa 

dX* 


=  3/=^=^. 


In  coordinates,  terms  on  the  right-hand  side  of  (100)  can  be  written 

F(X,  0  =  F“a ga  0  g“,  F(X,  0  =  ~F\ ~ga  ®  Ga . 


(101) 


(102) 


Basis  vectors  ga  and  their  reciprocals  g"  will  be  described  in  detail  later.  Both  F  and  F  are  second-order, 
two-point  tensor  fields  with  positive  determinants: 


detF  =  X-eahcea^FaMFbpFcmX  >  0,  detF  =  ^eaPxeABCFfAFbRFA.  >  0. 


1 

6l 


(103) 
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Figure  I.  Mappings  among  reference,  intermediate,  and  current  configurations  of  a  deformable  body. 


Permutation  symbols  are  ea^x  and  eapx;  relations  involving  these  implicitly  assume  n  =  3.  Inverting 

(100) — (102),  ' 

p-1  =  F-1F-1,  daXA  =  F-'am  =  F~lAaF~lama;  (104) 

t )  =  F~xamGa  <g>  g\  F-*(x,  0  =  F~l«a ga  <g>  g“.  (105) 


Furthermore,  from  the  definition  of  the  inverse, 


V 


-  1A 


=  8 


B  ’ 


i?a  77—  ltr  _  cfl 

mb  —  8b, 


pa  p~XA 
r  Jtr  .p 


—  pa  p~la  —  §a 


(106) 


with  Kronecker  delta  symbols  <57.  The  target  space  of  F  and  F  1  is  intermediate  configuration  B.  Figure  1 
illustrates  deformation  mappings  between  reference,  intermediate,  and  current  configurations  entering  (100). 

A  multiplicative  decomposition  such  as  (100)  can  be  used  to  represent  various  physical  behaviors.  Usually, 
F  is  associated  with  locally  recoverable  elastic  deformation,  such  that  each  material  element  in  configuration  B 
is  locally  unloaded  or  stress-free.  In  geometrically  nonlinear  crystal  mechanics  [2,  3,  9,  11,  14,  25,  33,  42,  43], 
F  is  associated  with  lattice  (i.e.  thermoelastic)  deformation,  and  F  is  associated  with  plastic  (i.e.  dislocation 
slip  enabled)  deformation.  In  multiplicative  descriptions  involving  explicit  thermal  deformation  [14,  15,  25], 
the  total  compatible  thermoelastic  deformation  F  is  split  into  locally  recoverable  elastic  deformation  F  and 
stress-free  thermal  deformation  F.  In  theories  of  twinning  [19,  20],  F  represents  stress-free  twinning  shear.  In 
theories  of  porous  or  damaged  media  [18],  F  represents  volumetric  expansion  associated  with  voids.  In  theories 
of  growth  in  biomaterials  [16,  25],  F  represents  stretching  and  mass  changes  associated  with  tissue  structure 
evolution. 

A  decomposition  somewhat  alternative  to  (100)  enables  description  of  materials  containing  initial  residual 
stresses  and  various  kinds  of  defects.  In  the  continuum  theory  of  inhomogeneous  bodies  of  Noll  [6],  a  decom- 
position  F  =  (Vy)K-1  is  proposed,  where  K  may  denote  a  generally  anholonomic  mapping  from  a  residually 
stressed  intermediate  configuration  to  an  unstressed,  but  generally  disconnected,  reference  configuration.  The 
elastic  deformation  associated  with  stresses  superimposed  on  the  intermediate  configuration  is  represented  by 
the  integrable  (i.e.  compatible)  deformation  gradient  Vy,  where  y  is  the  motion  associated  with  this  elastic 
deformation  gradient.  Mapping  F,  like  K,  need  not  be  integrable.  Wang  and  Truesdell  [23]  applied  similar  ideas 
towards  finding  universal  solutions  for  certain  problems  involving  pre-stressed  laminated  elastic  bodies. 


3.1.2.  Anholonomicity  conditions  and  partial  differentiation.  It  is  assumed  that  F(X,  t)  is  at  least  once  continuously  differ¬ 
entiable  with  respect  to  referential  coordinates  XA.  Single-valued  coordinates  xa(X,  t)  referred  to  intermediate 
configuration  B  that  are  at  least  two  times  differentiable  with  respect  to  XA  exist  if  and  only  if  the  following 
integrability  conditions  apply  [27]: 


«  =  dBF« 


d2xa 

dXAdXB 


d2xa 

dXBdXA ' 


(107) 


When  (107)  applies  globally  (i.e.  for  allX  €  B0,  which  presumably  is  associated  with  a  simply  connected  body) 
intermediate  configuration  B  is  said  to  be  holonomic,  and  then  like  the  reference  configuration  B0,  is  labeled  a 
Euclidean  ^-dimensional  space: 


\aF:b]  =  0  O  =  3 Axa. 


(108) 
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Otherwise,  single-valued  coordinates  xa(X,  t )  continuously  (two  times)  differentiable  with  respect  to  XA  do  not 
exist  in  intermediate  configuration  B,  which  is  then  termed  an  anholonomic  space: 

d[A Fj]  ^  0  <=>  xa(X,  t)  anholonomic.  (109) 

The  domain  of  F  must  be  simply  connected  to  ensure  that  the  left  equality  of  (108)  constitutes  sufficient  con¬ 
ditions  for  existence  of  a  uniform  covering  xa  of  that  domain,  in  which  case  both  x"  and  XA  are  single  valued 
for  all  A  e  B0.  When  (108)  (respectively,  (109))  applies  only  over  local  simply  connected  regions  of  B0,  then 
deformation  map  F  is  designated  as  holonomic  (respectively,  anholonomic)  only  in  those  regions.  Regardless 
of  which  of  (108)  or  (109)  applies,  partial  differentiation  with  respect  to  intermediate  coordinates  is  defined  as 
follows  [5,  27,  28,  31]: 

3«0)  =  ■  (HO) 


Arguments  in  (107)— (1 10)  can  be  repeated  for  F  *(x,  t),  which  is  assumed  continuously  differentiable  with 
respect  to  spatial  coordinates  xa.  Single- valued  coordinates  xa(x,  t)  referred  to  intermediate  configuration  B  that 
are  continuous  and  at  least  two  times  differentiable  with  respect  to  xa  exist  if  and  only  if 


‘daF-\  =  dbF~ll  O 


2~a 


3-x' 


2~o/ 


dx' 


dxadxb  dxbdxa 

When  (111)  applies  globally  on R,  intermediate  configuration B  is  holonomic: 

3[ aF-\  =  0  4*  F~\  =  da~xa. 

Otherwise, 

d[aF~lamb]  /  0  xa(x,  t )  anholonomic. 

Partial  differentiation  with  respect  to  intermediate  coordinates  also  obeys 

3«(0  =  U-)K. 

Verification  that  (1 10)  and  (1 14)  are  equivalent  is  straightforward  using  (85)  and  (101): 

3«(0  =  M-)F-]Aa  =  da(.)  3AxaF-]Aa  =  90(.) FyP~XAM  =  3 a(-)F“a. 
Second  partial  anholonomic  derivatives  obey  the  relations 


(111) 

(112) 

(113) 

(114) 

(115) 


9*[3/>Q]  =  UM-)~F-l%]F-lAa  =  3^3a(-)]  +  2dA(-)d[a~F-XApv  (116) 

Following  similar  arguments, 

3«[3/i(-)]  =  3/i[3a(-)]  +  23fl(-)9[«Afl/i].  (117) 

In  general,  second  partial  anholonomic  differentiation  is  not  symmetric.  Explicitly, 

3[a[3ffl(-)]  =  -  3/>[3«G)]}  =  =  U-)\aF%v  (118) 

Only  when  (108)  and  (112)  apply  such  that  B  is  holonomic  is  second  partial  anholonomic  differentiation  is 
symmetric: 

F~lAmf}  =  3 PXA  =»  d[aF-Umfi]  =  d[a dp]XA  =  0  =►  3„[3^(-)]  =  3/J[3a(-)],  (1 19) 

F%  =  dpxa  =»  3 [aF^}  =  d[adP]xa  =  0=>  3ff[3^(-)]  =  3/?[3«(  )]-  (120) 


Even  when  (108)  and  (1 12)  apply  globally  on  Bq  and  B,  respectively,  inverse  coordinate  functions  XA(x,  t)  and 
xa(x,  t)  may  only  be  available  locally  [6]. 
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3.1.3.  Anholonomic  basis  vectors  and  metric  tensors.  By  definition,  basis  vectors  and  their  reciprocals  in  intermediate 
configuration  B  obey 

=  (121) 

A  symmetric  metric  tensor  g  on  B  is  defined  in  components  in  terms  of  a  scalar  or  dot  product  as 

gap  =  ga  ■  gfi  =  gft  ■  ga  =  gfta  =  g(afi).  (122) 

The  dot  product  of  two  generic  contravariant  vectors  V  =  Vaga  and  W  =  Waga  is  computed  as 

v  •  w  =  vaga  ■  wp gp  =  r  wp( ga  ■  gp)  =  vagaP  w p  =  vawa  =  va  wa.  (123) 

As  indicated,  the  metric  tensor  can  be  used  to  lower  contravariant  indices  in  the  usual  manner: 

Va  =  VPgap,  ga  =  gafigft.  (124) 

Metric  g  is  assumed  to  be  positive  definite,  with  positive  determinant  g  over  any  volume  (i.e.  excluding  possible 
points,  lines,  or  surfaces  where  g  may  be  zero  or  undefined): 

g  =  det  g  =  det(g„^)  =  \eafixeS€(l,ga&  gp€  gx<j)  >  0.  (125) 

D 

The  inverse  g_1  with  components  g“^  on  B  obeys,  by  definition, 

gTp  =  r  .  ~gP  =  ~gP  .  r  =  g*  =  g“*gx/;  =  ^  (!26) 


The  inverse  metric  (126)  enables  the  dot  product  of  generic  covariant  vectors  a  =  aa  g“  and  ft  =  paga  on  B: 

a  ■  0  =  aa g“  •  Ppgr‘  =  aaPfj( g“  •  gp)  =  aagaPPp  =  aa =  aapa.  (127) 

Components  of  g_1  can  be  used  to  raise  indices  as  shown  above.  Permutation  tensors  in  configuration  B  are 
defined  as 

€aPx  =  -4  <?Px,  eapx  =  y/§  eapx.  (128) 

Vg 

Shifter  tensors  can  also  be  introduced  among  basis  vectors  in  intermediate  and  reference  or  spatial  configura¬ 
tions: 


g 


.a  A  ~a 


a  a  ~a  ^jaa 


Furthermore,  since 


>;  g“  =  (g°3gfl>>  g£  =  < g«,ga>; 

(129) 

>  &aA  —  So*  '  (-*A  —  ^Aan 

(130) 

>  Soca  —  *  §«  —  &aa  ? 

(131) 

K,  n!igj  =  £&  =  *}; 

(132) 

det  (gap)/  det  (Gab)  =  \/g/G; 

(133) 

/ det  (gah)/  det  (gap)  =  yg/I- 

(134) 

ga  =gaAGA=gaaga; 

(135) 

1  •  a  —  ara  a  Qa  —  gra  Qa 

(136) 

i  =  g^GA  O  g *  =  g“g^, 

(137) 
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it  follows  from  the  product  rule  of  determinants  that  det  (g£)  =  [det  (g“)][det  (gA )]  =  [y/g/g][y/g/G]  =  y/g/G 
in  agreement  with  (26). 

First  consider  the  case  when  (108)  applies,  such  that  B  can  be  regarded  as  a  Euclidean  space  and  a  position 
vector  x  can  be  assigned  to  any  point  x(X,  t)  e  B.  In  that  case,  basis  vectors  can  be  defined  in  the  usual  manner, 
similarly  to  (4): 

ga(x)  =  dax=^.  (138) 

dxa 

Metric  tensor  components  corresponding  to  (138)  are 

gap(x)  =  3«x  •  dpi.  (139) 

Any  time-independent  Euclidean  coordinate  system  (e.g.  fixed  curvilinear  or  Cartesian  coordinates)  can  be  used 
for  x"  in  this  case. 

Next  consider  the  case  when  (109)  applies,  such  that  single-valued  coordinates  xa(X,  t)  continuously  dif¬ 
ferentiable  with  respect  to  XA  do  not  exist  in  B.  In  that  case,  a  one-to-one  correspondence  between  a  material 
particle  X  and  an  intermediate  point  x  is  not  applicable.  Intermediate  basis  vectors  associated  with  a  given 
material  particle  should  then  be  assigned  to  X  rather  than  3c,  i.e.  ga(X).  Any  time-independent  coordinate  sys¬ 
tem  could  be  used  for  ga(X)  in  this  case  (time-dependent  convected  systems  are  considered  later).  However, 
it  may  prove  practical  to  select  identical,  possibly  curvilinear,  coordinate  systems  in  configurations  B0  and  B 
[14,36]: 

~ga(X)  =  SA  Ga  =  SaJaX  =  SiF^dpX.  (140) 

Metric  tensor  components  corresponding  to  (140)  are 

gap(X)  =  ga-gp  =  8aGa  •  8f Gb  =  SA8BGAB,  (141) 

where  GAB(X)  is  the  referential  metric  tensor  of  (10).  Mixed  Kronecker  delta  symbols  are  equivalent  to  shifter 
tensors  for  coincident  coordinate  systems  and  have  the  usual  meaning,  i.e. 

SA  =  IVA  =  a,  8a  =  OVA  ±  a.  (142) 

The  determinant  of  the  metric  is  simply 

g(X)  =  det(gff/j)  =  detfit)  det(3j)  det  (GAB)  =  det  (GAB)  =  G.  (143) 

Explicitly,  shifter  tensor  components  are 

gaA  =  <r ,  G^>  =  (8bGb,  Ga)  =  8°8ba  =  SaA,  (144) 

with  determinant  (133)  reducing  to 


det(g^)  =  1/  det(g“)  =  yf&jG  =  1.  (145) 

Now  consider  the  case  when  (113)  applies,  such  that  single-valued  coordinates  xa(x,  t)  continuously  differ¬ 
entiable  with  respect  to  xa  do  not  exist  in  B.  In  that  case,  a  one-to-one  correspondence  between  a  spatial  point 
x  and  an  intermediate  point  x  is  not  applicable.  Intermediate  basis  vectors  associated  with  a  given  spatial  point 
should  then  be  assigned  to  x  rather  than  x,  i.e.  g„(x).  Any  time-independent  coordinate  system  could  be  used  for 
go,(x)  in  this  case.  However,  a  pragmatic  choice  corresponds  to  identical  coordinate  systems  in  configurations  B 
and  B  [14]: 

&»(*)  =  S"ga  =  8aa  dax  =  8^F~^adpX.  (146) 

Choices  (140)  and  (146)  differ  but  are  not  contradictory;  in  particular,  basis  vectors  in  (140)  and  (146)  are 
related  by  8°ga[x(X,  /)]  =  8“gA (x,  A )SA  SB GB(X).  Metric  tensor  components  corresponding  to  (146)  are 

gat i(x)  =  ga  ■  gp  =  8aaga  ■  8bp gb  =  8aJbgab,  (147) 
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where  gab(x)  is  the  spatial  metric  tensor  of  (1 1).  Mixed  Kronecker  delta  symbols  satisfy 

=  Via  =  a,  8aa  =  OVa  /  a. 


(148) 


The  determinant  of  the  metric  is 

g(x)  =  det(  gafj)  =  dct((5")  dct(V^)  dct(gfl/,)  =  det(gflfe)  =  g.  (149) 

Explicitly,  shifter  tensor  components  are 


gao  =  (ga,ga)  =  (8bagb,ga)  =  8bA  =  K, 

with  determinant  of  (134)  reducing  to 

det(g“)  =  1/  det(g")  =  Vg/I  =  1. 


(150) 


(151) 


Finally,  the  simplest  choice  of  coordinate  system  for  configuration  B  is  a  Cartesian  system  with  constant 
basis  vectors  e(y : 


ga  —  gap  —  ^a  '  ^/3  —  ^of/3 ,  g  —  dct( 8ap )  —  1* 


(152) 


When  anholonomicity  conditions  (109)  or  (113)  apply,  such  a  Cartesian  frame  is  assigned  externally  to 
configuration  B,  since  in  that  case  B  is  not  a  Euclidean  space. 


3.1.4.  Convected  anholonomlc  connection  coefficients.  Arbitrary  connection  coefficients  in  reference  configuration 
Bo  can  be  mapped  to  coefficients  f ""  in  configuration  B  via  [14,  27] 


p..a  _  pa  £— 1 B  p-ic  ,  pa  3  p-lA  _  pa  £—15  £— 1C  _  £—15  p-\A  a  pa 

lPx~r-Ar  -r  .xlBC^rA°Pr  .x~r^r  .pr  ■X15C  r  .pr  .X°Br  -4’ 


^  ./T 


(153) 


where  (110)  applies.  Torsion  tensor  components  TRR:  of  (32)  in  configuration  Bo  map  correspondingly  as 


rhmmCl  T^CX  TJ—\B  TJ—\C  JimmA  jjCl  JJ—XB  Tf—XC  T''  m.A  p  ..O'  i  ..O'  ^..Q ' 

1  Px  -  -r  -X^C  -  -x1  [BC]  -  1  \fix\  KPx  -  1\fixV 


(154) 


where  components  of  the  anholonomic  object  k  are  defined  as  [27] 


j\..a  _  p—lA  £—15  c,  pa  _  j\..a 
KPX~P  -PP  .x0[aP -B]  ~  K[PX I' 


(155) 


Note  also  the  identities 


'■ap 


=  F 


r— 1 A  p'—XB 


\a~F % 


=  F 


—  \A 


■[“ 


F-w„dAF^  = 


p-lA  px  a,/?-1"  — 

r  .\ar  .\B\°Ar  .fll  — 


-15 


•PV 


■[<**  -|5| 


•P] 


FBd[a 


F 


-15 


■PV 


(156) 


When  F(X,  t)  is  an  integrable  function  of  XA  and  hence  (108)  applies,  then  =  0.  Assuming  now  that  F 

is  at  least  twice  differentiable  with  respect  to  reference  coordinates  XA,  Riemann-Christoffel  curvature  tensor 
components  R"R'A:n  of  (33)  in  configuration  B0  map  as  [14,  27] 


hmmma  _  pa  p—  XB  p— 1C  p—  XD  n...A  _  ^  p ..o'  ,  ^p 

K0X&  ~  * ^  -r  -X*  .8KBCD  ~  2d[P[  X]S  +  21 

Relations  (154)  and  (157)  are  consistent  with  (34)  and  (35)  [27]: 


p..e  -1-  _  J?mmM 


V[/jvx]r  =  vw(az]r  +  r^vs)  =  1-rj*s 

v[/svx]“5  =  Vw(3x]ati  -  f"“3a„)  =  --Rj*8 


vs-ffxvsva, 

Ota 


(157) 


(158) 

(159) 
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where  Va  and  as  are  respective  anholonomic  components  of  twice-differentiable  vector  and  covector  fields, 
and  covariant  differentiation  with  respect  to  anholonomic  coordinates  is  defined  formally  later  in  (179).  Skew 
components  of  second  partial  derivatives  in  (1 18)  can  be  written  as  follows  in  terms  of  the  anholonomic  object: 

3[«[3/i](-)]  =  =  dx(-)F~] (160) 

Relationships  analogous  to  (153)— (157)  can  be  used  to  map  connection  coefficients  and  torsion  and  curva¬ 
ture  tensors  referred  to  spatial  coordinates  to  intermediate  space  B  by  replacing  referential  indices  with  spatial 
indices  (e.g.  XA  ->  x“)  and  replacing  components  F\  with  F~Xo[a,  the  latter  assumed  twice  differentiable  with 
respect  to  xa . 

A  convected  coordinate  representation  of  F(X,  t)  can  be  used  to  verify  transformation  (153)  when  = 
P  corresponds  to  the  Levi-Civita  connection  of  (42)  for  Euclidean  space  on  Bo.  Convected  anholonomic 
basis  vectors  and  their  reciprocals  are  defined  as 


g'JX,  t )  =  F-lAa Ga,  g"*(X,  t)  =  F\ Ga- 

Cg'a^)  =  F~1AaF'B{GA,  Gb)  =  F-]AaF^8BA  =  F-'iFti  =  8%. 


The  deformation  map  of  the  second  of  (102)  can  then  be  expressed  as 

F(X,  t )  =  F\ ga  0  Ga  =  ga  0  F^Ga  =  ga  0  g'“  ■ 
Components  of  the  metric  tensor  corresponding  to  basis  vectors  in  (161)  are 


g  a-gp 


=  F 


—  1^  C'  T7~  m  r-  _  t?—  \A  rr—lB  r-  r*  _  77—  IA  /-■  77 

*J.4  •  r  [Ajb  —  r  ,ar  a^A  '  —  r  ,a(~rABX 


-IB 


(161) 

(162) 


(163) 

(164) 


Taking  the  partial  anholonomic  derivative  of  g'a  using  (110)  results  in 

3,sia  =  MF-LlGA)  =  Gr  fA  -  f-x%f-'ajbf'a)  g;  =  r~X 

in  agreement  with  (153).  The  torsion  tensor  of  (154)  vanishes  identically  in  this  case  from  (42): 

T’-.Q'  rfcx  rp—lB  t? —  1 C  G  ..A  p ..o'  .  ~..a  r\ 

bx  ~  b-Ab  -Pb  -XT  [BC]  ~  1  [fix]  +  KPX  ~ 
where  the  anholonomic  object  is  the  same  as  that  of  (155): 


~..ot  _  p- IA  p-lB  a  pa 
KPx  - r  -Pr  .x°iAr-By 


(165) 


(166) 


(167) 


Following  an  analogous  approach  with  F  1  (x,  t),  convected  anholonomic  basis  vectors  and  their  reciprocals 


are 


gjx,  t)  =  F“aga,  g ,a(x,t)  =  F  ‘“gfl. 

(g'a,  g'P)  =  KF~lf>.b( ga,  gh)  =  ~FaJ~Va  =  ~FaJ^.a  =  C 


(168) 

(169) 


The  deformation  map  of  the  first  of  (102)  can  be  expressed  as 

F(A,  t)  =  F“aga  cg>  g“  =  ga  0  ga. 

Components  of  the  metric  tensor  corresponding  to  basis  vectors  in  (168)  are 

ga  •  g^j  =  K ga  •  F% gh  =  KF^pga  ■  gb  =  F“agabFhmfi 
Taking  the  partial  anholonomic  derivative  of  g'a  using  (114)  results  in 

bZ'a  =  WKaga)  =  ba  ~  g' 


X  1  P<*  &X 


(170) 

(171) 


(172) 
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The  torsion  tensor  of  (154)  vanishes  identically  in  this  case  from  (53): 

.  . . _  +  <  =  0, 

where  the  corresponding  anholonomic  object  is 


' T’..Qf  77—  la  jjb  -r'C  g  ..a  p. 

mP  mX  r  M  “  [0X] 


-..O'  _  770  776  3  p—la 

Kfix  -  r.fir.ya[«r  .by 


(173) 


(174) 


■il  of  (165)  and  T"* 

decomposition  (101),  it  can  be  shown  that  their  skew  covariant  components  are  equal  [14,  p.  651]: 


Connection  coefficients  f  of  (165)  and  T’j*  of  (172)  are  generally  different.  However,  following  from 


r“x  —  —  F°  Fa  r)r,F  *  —  —F~ 
1  [/3a]  —  r.pr.a°[br  .a]  —  r 


■P 


F~ 


AbF^  —  r[/ja]. 


Thus,  anholonomic  objects  of  each  connection  defined  in  (167)  and  (174)  are  equal: 


-■■x 


—  _r"*  —  _r"x  — 

—  1  r«/vi  —  1  r«/vi  —  K- pa  • 


“-Pet  ~  1  [pa]  ~  1  [fia] 

Computing  skew  partial  derivatives  of  convected  basis  vectors  as 


=  T[afs\£x  =  =  Kpl  g'x, 


the  following  local  integrability  conditions  are  equivalent: 


9[ag«  =  0  ^  3[ag'p]  =  0  =  0  O  F\  =  dAxa  &  F-\  =  3a. 


=  OnX  . 


(175) 


(176) 


(177) 


(178) 


3.2.  Anholonomic  covariant  derivatives 

Differentiation  with  respect  to  general  anholonomic  coordinates  is  developed.  Covariant  derivatives  and  cor¬ 
responding  connection  coefficients  are  defined.  Various  choices  of  anholonomic  connection  coefficients  corre¬ 
sponding  to  different  basis  vectors  in  the  intermediate  configuration  are  examined.  Total  covariant  derivatives 
of  two-point  (and  three-point)  tensors  with  one  or  more  indices  referred  to  anholonomic  space  are  defined. 
Divergence,  curl,  and  Laplacian  operations  and  corresponding  identities  are  presented. 


3.2.1.  Differentiation.  Covariant  differentiation  with  respect  to  anholonomic  coordinates  is  defined  similarly  to 
(31).  Let  A  be  a  vector  or  higher-order  tensor  field  with  components  AJy^ .  The  covariant  derivative  of  A  is 
computed  as 


U  AU-<t>  _  Q  AOt—<t>  I  T'..ot  AP--4  _i_  .  .  .  _|_  p..0  xa...p 
VVAy...p  —  °VA  y...p  T  1  vp  Ay...p  +  +  1  vpAy...p 


p..pAa,,.(/>  _  _  _  _  _  p..p^a,..</> 

vy  p...p  vp  y...p- 


(179) 


Partial  differentiation  obeys  (115).  Connection  coefficients  referred  to  intermediate  space  B,  written  as  T"“, 
in  general  consist  of  up  to  n 3  entries,  where  n  is  the  dimensionality  of  Euclidean  spaces  B o  and  B.  Partial 
derivatives  of  intermediate  basis  vectors  and  their  reciprocals  obey  (restricting  attention  to  metric  connections), 
by  definition, 

9/3ga  =  r^gx,  9/ig“  =  -r^g*.  (180) 

Therefore,  covariant  derivatives  of  intermediate  basis  vectors  vanish  identically: 

V/3 ga  =  9/iga  -  r^gx  =  0,  V^g“  =  dp g“  +  r-g*  =  0.  (181) 

For  example,  applying  (179)— (180)  to  generic  vector  field  A  =  A01  g„  gives  the  gradient  of  A: 

dpA  8  g^  =  9 p(Aaga)  8  gp  =  (3pAa+Ax r-^)ga  8)  g^  =  VpAa g„  8  g*.  (182) 

By  definition,  the  gradient  of  a  scalar  field  is  the  same  as  its  partial  derivative,  e.g. 

VaA  =  daA  =  F~lidBA  =  FbadbA.  (183) 
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Recall  from  (176)  that  the  anholonomic  object,  denoted  generically  by  .7*  =  =  .7*  ,  obeys 


7— 1 B  77 —  1A 


KPa  —F~iDaF~ 


1  77— M  c\  77’ X  77 X  C\  I?—IA  T-'b  77 £7  O  77  AX  77  AX  O  77#  ,"X 

,p^  ^  m]  ~  *  .fif  M0[br  .a]  —  ^  .a°[pr  .a]  ~  K[Pa]- 


<b  rya 


7-1/ 


(184) 


The  intermediate  torsion,  denoted  generically  by  T"*,  is  defined  using  the  final  two  equalities  in  (154): 


T"a  —  r-«  i  i,--01 

LPx  -  1  IPx]  +  KPx 


_  rT~'mmO( 

~  1lPxY 


(185) 


Here,  the  first  two  equalities  in  (154)  are  not  required  to  hold  for  this  generic  definition  of  the  torsion  of  an 
anholonomic  space  (i.e.  the  torsion  in  (185)  need  not  map  between  configurations  as  a  true  tensor).  Similarly, 
the  intermediate  curvature,  denoted  generically  by  Rp“s,  is  defined  using  only  the  final  two  equalities  in  (157): 


R7xs  =  2«  +  +  2.-17“  = 


(186) 


Here,  the  first  two  equalities  in  (157)  are  not  required  to  hold  for  this  generic  definition  of  the  curvature  of  an 
anholonomic  space  (i.e.  the  curvature  in  (186)  need  not  map  between  configurations  as  a  true  tensor). 

With  definitions  (184),  (185),  and  (186)  now  given,  skew  gradients  can  be  obtained.  Specifically,  from  (160), 


9[o'[9/3](-)]  —  Kap^x^y 


(187) 


For  the  twice-differentiable  scalar  field  A, 


V[a(Vp]A)  —  (![ff(3/j]4)  r{*p]dxA  -  (Kap  +  r^)3xv4  —  TapdxA. 


(188) 


It  can  be  shown  that  (158)  and  (159)  also  hold,  i.e. 

vpvx]r*  =  1-RmV  -  Tfyvsv°, 


lPxS' 


Px 


v[/svx]a<5  =  --Rj“saa  -  Tj£Vaas 


(189) 


where  Va  and  ag  denote  components  of  twice-differentiable  vector  and  covector  fields,  respectively. 


3.2.2.  Anholonomic  connection  coefficients.  Particular  choices  of  coefficients  T^“  entering  (179)  and  (180)  are 
discussed  in  what  follows. 

First  consider  the  case  when  (108)  applies,  such  that  B  can  be  regarded  as  a  Euclidean  space  and  a  position 
vector  x  can  be  assigned  to  any  point  x(X,  t )  e  B.  In  that  case,  basis  vectors  and  metric  tensor  components  can 
be  defined  as  in  (138)  and  (139),  and  any  time-independent  Euclidean  coordinate  system  (e.g.  fixed  curvilinear 
or  Cartesian  coordinates)  can  be  used  for  xa.  Coefficients  of  the  Levi-Civita  connection  on  B  are,  analogously 
to  (42), 

r£(*)  =  f  fx  =  \sa\dp~gx «  +  3 xgps  ~  3 sgpx)  =  f  (190) 

Since  xa  are  holonomic  coordinates  in  Euclidean  space,  the  torsion,  anholonomic  object,  and  curvature 
associated  with  coefficients  (190)  all  vanish  identically.  Skew  gradients  of  basis  vectors  also  vanish  [28]: 
3[ag£]  =  3[a3,8]X  =  0. 

Next  consider  the  case  when  (109)  applies,  such  that  single-valued  coordinates  xa(X,  t)  continuously  differ¬ 
entiable  with  respect  to  XA  do  not  exist  in  B.  When,  following  (140),  identical  coordinate  systems  are  used  in 
configurations  B 0  and  B: 

ga(X)  =  8AGA,  (191) 

rw  =  =  8ac8pDGCD8AGA  =  8“GcaGa  =  8°  GA,  (192) 

gap{X)  =  8a8*Gab,  ~ga\X)  =  8aA8pBGAB.  (193) 
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Taking  the  partial  derivative  of  (191)  and  applying  (45)  and  (110)  gives 

dpg*  =  8%GA  =  S^F~yB Ga  =  8AaF~w£  fAGc  =  £s*~gx  =  r£gx.  (194) 

Similarly,  taking  the  partial  derivative  of  (192)  and  applying  (46)  and  (110)  gives 

3^g“  =  SaAdpGA  =  S“F~l *pdBGA  =  -8aAF~XB/T  icGC  =  SaAF~lB/r  fc8cx~gx  =  -T£g*.  (195) 

Comparing  (180),  (194),  and  (195),  connection  coefficients  consistent  with  (140)  are  time  dependent: 

rfx(X,  t)  =  858^  ^c(X)~F-wp(X,  t).  (196) 

Since  Levi-Civita  connection  coefficients  p  are  symmetric,  it  follows  that  r defined  in  (196)  obeys 

■p.-a  _  ooc  oCG  mmA  r?—  IB  _  oa  oCG  ••A  rp—  \B  _  ooc  oBG  m.A  77—  1 C 

1  —  °A°xY  BC*^  .^—^a^xTcb^  .p  —  °A°xTCB*  ■£’  V W!) 


Covariant  indices  of  T'F  are  generally  not  symmetric;  the  left  covariant  component  corresponding  to 
differentiation  by  Vp(-)  in  (179)  correlates  with  F~ 1  [p .  The  torsion  of  (196)  is  defined  as  in  (185): 


7 i  z-"a 

LPx  —  1  VPx]  ^  KPx 


_  ?0(G  mmA  77 

—  °A  r  BC r 


-IB  oC 


<Ni  +  F 


—  IA 


■lPux] 


■IP 


F-W'X]dAF[ 


.B’ 


(198) 


where  i<Px  is  the  anholonomic  object  associated  with  F  of  (167).  When  the  anholonomic  object  k"“  is  non-zero, 
skew  partial  derivatives  as  in  (187)  are  generally  non-zero.  When  the  torsion  Tpx  is  non-zero,  skew  covariant 
derivatives  of  a  scalar  field  as  in  (188)  are  generally  non-zero.  The  Riemann-Christoffel  curvature  associated 
with  (196)  is  defined  as  in  (186): 


=  29«  +  2T- ,T-  +  2ic px T"“  =  -25 


G  .Jt 


a?Dp—lB  &-1C  o  L 
°*r  .[Pr  .x]°BT  CD’ 


A  8  ■ 


(199) 


where  the  vanishing  of  the  curvature  tensor  of  Euclidean  reference  space  (43)  has  been  used.  When  curvature 
R‘PxS  of  (199)  and  torsion  77"  of  (198)  do  not  vanish,  skew  covariant  derivatives  of  vector  and  covector  fields  are 
generally  non-zero,  as  in  (189).  It  is  noted  that  torsion  and  curvature  defined  in  this  way  can  each  be  non-zero 
even  when  configuration  B  is  holonomic,  i.e.  even  when  E"  —  Ti.r"  and  k"“  =  0.  Using  (49),  the  (negative) 
covariant  derivative  of  the  metric  tensor  is 


Maftx  —  ^agPx  ~  dagpx  +  1 ^apSSx  +  1", 


ax 


gps  =  -8Bp8cxF-liGBC;A  =  o. 


(200) 


Therefore,  anholonomic  covariant  differentiation  commutes  with  lowering  indices  via  the  metric  gap  = 
8Aa8BGAB. 

Now  consider  the  case  when  (113)  applies,  such  that  single-valued  coordinates  xa(x,  t)  continuously  dif¬ 
ferentiable  with  respect  to  xa  do  not  exist  in  B.  Following  (146),  identical  coordinate  systems  are  used  in 
configurations  B  and  B: 

&»(*)  =  8aa  ga,  (201) 

r  (X)  =  ~ga%  =  8ac8pdgcds;ga  =  8?g™ga  =  8“  g“,  (202) 

gap  (x)  =  8aa8bgab,  ~g^(x)  =  8aXg“b.  (203) 


Taking  the  partial  derivative  of  (201)  and  applying  (56)  and  (1 14)  gives 

dpga  =  8aJp ga  =  8aaFbp dbga  =  8"aFb/T  tag c  =  8aaF\\  £8* ~gx  =  T^gx. 
Similarly,  taking  the  partial  derivative  of  (202)  and  applying  (57)  and  (114)  gives 


dpga  =  8°  dpg“  =  8«aFbdbg“  =  —8aaFbJL  £gc  =  -8aaFb*  £8'g*  = 


1  Px%  ■ 


(204) 


(205) 
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Comparing  (180),  (204),  and  (205),  connection  coefficients  consistent  with  (146)  are  time  dependent: 

r?x(x,  t )  =  t(x)Fh'P(x,  t).  (206) 

Since  the  Levi-Civita  connection  p  y“  is  symmetric  in  covariant  indices,  it  follows  that  T of  (206)  obeys 


_  Oi 

r4x  -  5 


a  oc  g 


■:“Fba  =  8“8 


xr bc  -p 


XT  cb‘  .fi 


'Fa  =  8‘ 


tv bg  ••a tjc 

°xTcbr4 


(207) 


Covariant  components  of  V"“  are  generally  not  symmetric;  here  the  left  covariant  component  corresponding  to 
differentiation  by  V^(-)  in  (179)  correlates  with  F\.  The  torsion  of  (206)  is  defined  as  in  (185): 


j..a  _  p..ff  ,  -..a 

lPx  ~  1  [fix]  +  KPx 


_  sag  ..a-r'b  pc  ,  r’H  77 

—  0„  ^  h„P  r«0„i  +  r 


of bc  44  x] 


IP  -X] 


°ar  .b’ 


(208) 


where  /?"')  is  the  anholonomic  object  associated  with  F  of  (174).  When  anholonomic  object  is  non-zero, 
skew  partial  derivatives  as  in  (187)  are  generally  non-zero.  When  torsion  T""  is  non-zero,  skew  covariant  deriva¬ 
tives  of  a  scalar  field  as  in  (188)  are  generally  non-zero.  The  Riemann-Christoffel  curvature  of  (206)  is  defined 
as  in  (186): 


R 


PX‘ 5 


=  23[/jT"“+2r-r-  +2^r;-“ 


=  -28' 


g  .. a 

P  cd  ’ 


(209) 


where  the  vanishing  of  the  curvature  tensor  of  Euclidean  current  space  (54)  has  been  used.  When  the  curvature 
R"F  of  (209)  and  torsion  T"“  of  (208)  do  not  vanish,  skew  covariant  derivatives  of  vector  and  covector  fields 
are  generally  non-zero,  as  in  (189).  Note  that  the  curvature  and  torsion  defined  in  this  manner  need  not  vanish 
even  if  configuration  B  is  holonomic,  i.e.  even  if  F~lama  =  daxa  and  anholonomic  object  kmF  =  0.  Connection 
coefficients  and  curvature  for  this  case  are  worked  out  explicitly  for  the  choice  of  cylindrical  coordinates  in  the 
Appendix.  The  (negative)  covariant  derivative  of  the  metric  tensor  is 

Mapx  =  -VagPx  =  -3  agpx  +  rfpgSx  +  r  aXgps  =  -&h^cxFc‘agbc-,a  =  0,  (210) 

where  (62)  has  been  used.  Thus,  covariant  differentiation  commutes  with  lowering  indices  via  gap  =  Kfygab- 
Finally  consider  the  simplest  case  whereby  Cartesian  bases  are  used  for  B,  as  in  (152): 

go  =  ea,  g"  =  e“,  gap  =  ea  ■?£  =  8aP,  gafi  =  e“  •  =  8afi.  (211) 

In  this  case,  partial  derivatives  of  basis  vectors  vanish  identically, 

dpga  =  dp  e„  =  F-lB'pdBea  =  Fbpdbea  =  r"*gx  =  0,  (212) 

as  do  partial  derivatives  of  their  reciprocals, 

dpga  =  df,ea  =  F-'%dBea  =  Fbp  dbea  =  -T^g*  =  0.  (213) 


It  follows  that  connection  coefficients  on  B  must  also  vanish  in  this  case: 


T"“  —  0 
1  Px  ~  u- 

Note  however  that  the  anholonomic  object  need  not  vanish;  hence,  the  torsion  of  (185)  is 


T"a  _  ..a 

LPx  ~  KPx- 


(214) 


(215) 


Curvature  tensor  (186)  vanishes  identically  since  T^“  =  0.  The  covariant  derivative  of  the  metric  also  vanishes: 


Mapx  —  V agpx  —  da8px  —  F  ,adA8px  —  Fmada&px  —  0- 


(216) 
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3.2.3.  Total  covariant  derivatives.  Covariant  differentiation  of  two-point  tensor  fields  with  one  or  more  components 
referred  to  anholonomic  coordinates  (i.e.  one  or  more  indices  referred  to  configuration  B)  is  defined  following 
arguments  similar  to  those  of  Section  2.3.2.  First  consider  a  generic  two-point  tensor  of  the  form  A(X,  t),  with 
components  AG  hM-  The  total  covariant  derivative  of  A  is  calculated  as 


I  p..a  AP-4  A...F  , 

“r  1  vpAy...p  G...M  T 

1  T-4  401-P  A...F 
“T  1  vpAy...p  G...M 

—  4a-4  4...F  _  _  _  _ 

1  vy^1p.„ii  G...M 

—  Y-P 4<x~4  A...F 

1  vp^y...p  G...M 

i  77— 1 NG  -A  ja...tp  R...F 
‘  r  .vY  NR^y-F  G...M 

1  1  fp—ING  m»F  Ad. ..(f)  A...R 

~r  '  '  '  ~r  r  ,vp  NRr*~y...p  G...M 

77— 1  NG  .. R  Aa...(j)  A...F 
r  .vY  NGAy...p  R...M 

rp—lN  G  ,.R  Ad.. .(f)  A...F 
r  .yp  NM^y...p  G...R 

1  G  mmA  Ad. ..(f)  R...F  1 

'  p  NR^y.-.p  G...M  '  ’  ’ 

1  G  ..F  Ad. ..(f)  A...R 
p  NR^y.-.p  G...M 

G  m.R  aci ...(J)  A...F 
r  NG^y.-.p  R...M 

G  wR  Ad. ..(f)  A...F~\rp—\N 

P  NM^y.-.p  G...RV  .y 

,  ap-4  a-f  1 

r  1  vpAy...p  G...M  T 

1  4<*-P  A-F 

“T  1  vpAy...p  G...M 

—  F"^ 4a-4  A...F  _  .  .  . 
1  vy^1  p...p  G...M 

/  Aa--4  A...F  \  Ty—IN 
~  (Ay...p  G...m);Nf 

_  F-P 401-4  A...F 

1  vp^y.^p  G...M 

,  4P..4  A...F  , 

r  1  vpAy...p  G...M  ^ 

1  401-P  A...F 

“T  1  vpAy...p  G...M 

p..p  Ad. ..(f)  A...F 

1  vy^p...p  G...M 

(  Ad.. .(f)  A...F  \  fp—\N 

~  \Ay...pG...MKN r  .y 

_  r-p  4a-4  a...f 

1  vp^y^.p  G...M 

(217) 


Next  consider  a  generic  two-point  tensor  of  the  form  A  (x,t),  with  components  The  total  covariant 

derivative  of  A  is  computed  analogously  to  (217),  replacing  referential  coordinates  XN  with  spatial  coordinates 
xn  and  F~  lNmV  with  F"v : 

(4<x~4  a--f  \  _  a  (  Aa-4  a-f  \ 

v/1y...p  u^\/±y ...p  g...m) 


1  AP-4  a-f  ,  . 

'  1  vp^y ...p  g...m  ' 

. .  _i_  F"1^ Aa  -p  a-f 

'  1  vp  ^y  ...p  g...m 

_  Y"PAa-4  a-f  _  . 
vy  p...p  g...m 

.  .  _  F-P  Aa-4  a-f 

vp  y...pg...m 

1  pn  S  ..a  Ad-4  r—f 
'  1  .yp  nr^y.-.p  g...m 

i  .  .  .  1  pn  g  ../  Aa-4  a...r 
'  '  1  .yp  nr^y...p  g...m 

_  p'n  g  ..r  Ad. ..(f)  a...f 

.yp  ng;y...p  r...m 

—  ...  —  Fng  Aa  -4  a-f 

.yp  nmrLy...p  g...r 

i) 

1  s  -a  Ad. .4  r...f  . 

'  p  nr'rLy...pg...m  ■ 

.  .  .  1  g  ../  Ad. . .(f)  a...r 

'  p  nr^y ...p  g...m 

_  g  ..r  Ad -4  a-f  _  . 

P  ngp^y.-.p  r...m 

..—g  ..r  Aa-4  a...gtpn 

P  nm^y...p  g...r V  »v 

1  AP-4  a..f  ,  _ 

'  1  vp^y ...p  g...m  ' 

. .  _i_  Y"^ Aa  "p  a  "f 

'  1  vp  ...p  g...m 

—  r-p  a a-f  —  ■ 

vy  p...p  g...m 

.  .  —  F-P Aa-4  a-f 
vp  y...p  g...m 
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.  (A<x  -<l>  a  -f  V  Z7« 

.v 

i  AP  ~<t>  a-f  i  . 
'  vp  y .../a  g...m  ' 

—  F"p Aa"  ^  a"f  —  • 

1  vy^p.../!  g...m 

.  (  ACt.,4  a...f  \  fyn 
G^y-P  g  ~m>-nr  .V 


.  .  Y-^AU  -P  a-f 
'  vp  y...pig...m 

.  .  _  P"P  Aa-<t>  a-f 
1  vpr^y  ...p  g...m 


For  example,  letting  A(X,  t)  =  F(X,  t)  =  F\ ga  <g>  GA,  the  material  gradient  is  computed  using  (217)  as 
yF  =  dBF  0  Gs  =  (a*F“  +  T  F*  -  f  ^.c)  g«  ®  Gs  =  F"  :Bga  ®GA  ®  GB.  (2 1 9) 
Note  that  unlike  yF  in  the  first  expression  of  (89),  skew  covariant  components  of  yF  do  not  necessarily  vanish: 

nA:B]  =  +  t-f^f*,  -  G  $a]f:c  =  a  [bf^  +  r-  ^f*  .  (220) 

From  (220),  the  material  gradient  yF  is  generally  symmetric  in  covariant  indices  only  when  both  (108)  holds 
and  symmetric  connection  coefficients  V =  r(”7,  ]  are  prescribed  on  B.  As  an  example  of  (218),  letting 
A(x,  t)  =  F_1(x,  t)  =  F~XaM ga  <8>  ga,  the  spatial  gradient  is  computed  as 


i  r..«c 

+  1  Pxr 


—  lx  Z7— 1^  S  »■£  77 —  \OL 


^(F-1)  =  a/,F-1®g  h=(dh~F-\ 

Skew  covariant  components  of  y(F_1)  do  not  necessarily  vanish: 


F  ‘i  -  f  bCaF  '“■)  ®  ga  ®  g6  =  ^  ®  g°  ®  g6-  (221) 


F 


—  \a 


■[«:*] 


=  +  T"°lF~xp<uF~lx 


g  ..c 


,a]  -  r  =  9i*F 


■«]  ^  £x 
- g  xf'-l  1 


■[F 


'—\a 
■a] 


,  p..a  p-v  77 

1  [/3X]^  -'F 


—  1/3  p-lx 


.// 


(222) 


From  (222),  the  spatial  gradient  y  (F  )  is  generally  symmetric  in  covariant  indices  only  when  both  (112)  holds 
and  symmetric  connection  coefficients  =  r"(V  are  prescribed  on  B. 

Finally,  consider  a  tensor  field  of  order  three  or  higher,  with  components  referred  to  all  three  configurations 
5,  B,  and  B0,  written  A\ G/FM.  By  extension,  the  total  covariant  derivative  of  field  A  i 


is 


(  Aa-<t>  a-f  A  -F  V  _  a  (  AOt...(j)  a..f  A...F  \ 
y-™-y...Li  g...m  G...M)'V  uv\JFLy...Li  g...m  G...M ) 


+ 

p-“  AP- 
1  i >pAy. 

.4  a. 
-P  g- 

../  A...F 
..m  G...M 

+  • 

■  ■  +  r"pAy: 

■P  a-f 
■P  g-111 

A...F 

G...M 

Y"PA 

a.. 

.4  a. 

../  A-F 

- Y"PAa- 

f  a-f 

A...F 

vy 

P- 

■P  g- 

..m  G...M 

vp^y- 

p  g...m 

G...M 

+ 

F"g 

r  .up 

..a 

nr 

Aa- 

Ay.. 

4  r..f  A. 
4 1  g...m  G 

..F 

...M 

+  ---+nf 

■■/  Aa 
1  nr  ^y 

„4  A. 
■Pg- 

..r  A. 
,.m  G 

...F 

...M 

- 

fyng 

^.yp 

..r  acl... 
ng^-y... 

/p  a-f  A.. 
,/x  r...m  G. 

,.F 

..M 

—  ...  —  png 
r  .vy 

..r  acl. 
nm  y . 

.4  a.. 

...p  g.. 

f  A.. 
.r  G. 

.F 

..M 

(223) 


,  p-\NG  ..F  Aot. ,.<p  a..f  A...R 
m  r  .up  NRAy...p  g...m  G...M 

_  f-ING  ..R  au...<P  a-f  A...F 
r  .1/p  NMAy...p  g...m  G...R 


i  fy—lNG  ..A  AC/.. 4  a...f  R...F  . 

+  r  .1 >Y  NRAy-\Lg.'..mG-M~F 

_  f-ING  ..R  a<*~4  a-f  A...F 
r  .i/P  NGAy...p  g—m  R...M 

—  /  ac/. -4  a...f  A...F  \ 

g...m  G...M' -nr  . v 

_  /  Aa...d>  a...f  A...F  \  ip— IN 

x*y ...pi  g...m  G...m)'N^  .v* 

Consider  a  differential  line  element  c/X  in  the  reference  configuration.  Such  an  element  is  mapped  to  its 
representation  in  the  intermediate  configuration  dx  analogously  to  (86)  via  the  Taylor  series  [14] 


dxa(X)  =  (F“ ) 


xdXA  +  ^(~f:a.m) 


xdXA  dXB  +  ^(FZ,,bc) 


x 


dXA  dXB  dXc  + 


(224) 
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where  components  of  the  total  covariant  derivative  of  F  are  given  in  (219),  and  components  of  F“ rBC  can  be 
obtained  through  iteration  of  (217).  Similarly,  consider  a  differential  line  element  dx  in  the  current  configuration. 
Such  an  element  is  mapped  to  its  representation  in  the  intermediate  configuration  dx  via 


dxa(x)  =  (F~laJ 


dxa  +  ^(F-\,h) 


dxadxh  +  ^(F-\,hc) 


dxa  dxb  dxc  + 


(225) 


where  components  of  the  total  covariant  derivative  of  F-1  are  given  in  (221),  and  components  of  F~lama.bc  can  be 
obtained  through  iteration  of  (218).  To  first  order  in  dX  and  dx,  and  using  (90)  and  (100),  a  standard  assumption 
is 


dx  =  PdX  =  F“‘F  dX  =  F-1  dx. 


(226) 


3.2.4.  Divergence,  curl,  and  Laplacian.  Definitions  that  follow  correspond  to  anholonomic  space  B,  for  which  (179) — 
(183)  apply  for  covariant  and  partial  differentiation.  The  divergence  of  a  contravariant  vector  field  V  =  Vaga 
is 

(V,  V)  =  tr  (VV)  =  (3„v,  r>  =  V„r  =  da  r  +  r$  vp.  (227) 

Note  that  the  analog  of  the  final  equality  in  (66)  does  not  necessarily  apply  here,  since  (52)  may  not  hold  for 
connection  coefficients  on  B  (i.e.  T£“  are  not  necessarily  Levi-Civita  connection  coefficients).  The  vector  cross 
product  x  obeys,  for  two  vectors  V  and  W  and  two  covectors  a  and  y3, 

V  x  W  =  eapx  VpWxga,  a  x  p  =  euPxappx ga.  (228) 

The  curl  of  a  covariant  vector  field  a  =  aa  ga  is  then  defined  as 

V  x  a  =  g“  x  da(apgp)  =  ga  x  gpVaoip  =  eaPxVaapgx  =  eaPxVpax ga.  (229) 

A  relationship  like  the  final  equality  in  (68)  does  not  necessarily  hold  here,  since  coefficients  T^“  are  possibly 
non- symmetric.  The  Laplacian  of  a  scalar  field /  is,  from  the  symmetry  of  inverse  metric  gaP, 

v2/'  =  ~gapyavpf  =  ~gapVPVJ  =  ~gaPVp(dJ)  =  g^Adpf).  (230) 


An  identity  such  as  the  third  equality  in  (69)  does  not  necessarily  apply  here.  Unlike  (70),  the  divergence  of  the 
curl  of  a  (co)vector  field  does  not  necessarily  vanish  in  an  anholonomic  space: 


(V,  V  x  a)  =  Va(eapxVpax)  =  VJ(1  /y/§)VPxVpax 


(231) 


Identity  (189)  has  been  used.  Similarly,  from  (188),  the  curl  of  the  gradient  of  a  scalar  field  may  be  non-zero  in 
anholonomic  coordinates: 

V  x  V/  =  eaPx  (Vp  Vxf)~ga  =  eapx  (T%  3 Yf)~ga .  (232) 

Identities  (231)  and  (232)  are  now  examined  for  the  four  particular  sets  of  connection  coefficients  T^“ 
defined  in  Section  3.2.2. 

When  configuration  B  is  a  Euclidean  space  and  (190)  defines  the  connection  coefficients,  then  the  right-hand 
sides  of  (231)  and  (232)  are  identically  zero. 

When  (196)  defines  the  connection  coefficients,  i.e.  when  identical  (possibly  curvilinear)  coordinate  systems 
are  used  in  reference  and  intermediate  configurations,  then  (231)  and  (232)  become 


1 


(V,  V  x  a)  =  —e^i^F-^F-^d^  £)as  -  (8*D° 
x  F  mS(d ootx  —  S€e8^ p  ’^cae)]. 


&+F~ 


■[“ 


F~ 


■P] 


9 aF*j,) 


(233) 


1 


V  x  V/  =  -—eaPx[(8^^cF~^a8 


lPux\ 


+  WFaGF. 


(234) 
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Definitions  (198)  and  (199)  have  been  used  for  torsion  and  curvature.  Recall  that  neither  need  be  zero  even  if  F 
is  integrable.  Vanishing  of  the  covariant  derivative  of  g  =  det(gap)  follows  from  (200). 

Analogously,  when  (206)  defines  the  connection  coefficients,  i.e.  when  identical  (possibly  curvilinear) 
coordinate  systems  are  used  in  current  and  intermediate  configurations,  then  (231)  and  (232)  become 


(V, 


1 

V  x  a)  =  — e' 

Vs 


vPx 


~  (Kr  Wtn  +fVf.pV°f  J 


(235) 


x  F.s(ddax 


8€8cg 


x  r dc 


£«<)], 


xV/=-l^[(fr 


a  Y  be 1  .[Puxi 


(236) 


Definitions  (208)  and  (209)  have  been  used  for  torsion  and  curvature.  Recall  that  neither  need  be  zero  even  if 
F-1  is  integrable.  Vanishing  of  the  covariant  derivative  of  g  =  det(gap)  follows  from  (210). 

Finally,  consider  the  case  when  (214)  defines  the  connection  coefficients,  i.e.  when  a  Cartesian  system  is 
used  in  the  intermediate  configuration,  such  that  Y'F  =  0  identically.  In  that  case,  (231)  and  (232)  become, 
respectively, 


(V,  V  x  a)  =  Va(€^Vflax)  =  ea^Va(V^ax)  =  e^xd[a(d^ax)  =  ea^Kfa dbax, 
V  x  V/  =  caf>x(Vp Vx/)g„  =  e“^[3^(3x]/)]ea  =  e“^[^3/]ea. 


(237) 

(238) 


Thus,  even  if  intermediate  connection  coefficients  vanish,  the  divergence  of  the  curl  and  the  curl  of  the  gradient 
need  not  be  identically  zero  in  anholonomic  coordinates. 


3.3.  Anholonomic  Jacobian  determinants 

3.3.1.  Jacobian  determinant  of  F.  Jacobian  determinant  J  provides  the  relationship  between  differential  vol¬ 
ume  elements  in  reference  and  intermediate  configurations.  Differential  volume  element  in  the  intermediate 
configuration  ( dV )  is 

J  dV  =  dV  =  VldV  dx2  dx 3  c  B.  (239) 

The  Jacobian  determinant  J\¥(X,  t),g,  G(X)]  is  [14] 


J  =  \eABCeafix  FJVf^FV  =  l~WJG  eABCeaPxF^BFjc  =  VIJg  detF  =  / 


det(F“).  (240) 


From  similar  arguments,  inverse  Jacobian  determinant  J  1  [F  l(x,  t),  G(X),  g]  =  1/J  is 

J- 1  =  yaPx€ABC  ~F~XA.aF~W.p~F~XC.x  =  JWs  detF"1  =  yj det (Gab)/ det (gap)  det (F~XAM).  (241) 


When  the  mapping  is  restricted  to  rigid  translation  (or  to  no  motion  at  all),  then  Ff  t  =  g"  [the  shifter  of  (129)], 
F  =  g“ ga  ®  GA,  and  J  =  v/g/G’  det(g“ )  =  1  follows  from  (133).  From  identity  (51),  it  follows  that 


3J 


3  F‘ 


=  JF 


-1 A 


1-1 A  ~J  V -A' 


^4 


3  F 


(242) 


Taking  the  divergence  of  the  first  of  (242)  results  in  [14] 


(Jf~1AJ:A  =  \tapx*ABC  (f^bfV.a+K.af!c)  =  ^ABC~f^f\b-.cy 


(243) 
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Recall  from  (220)  that  material  gradient  yF  is  generally  symmetric  in  covariant  indices  only  when  both  (108) 

holds  and  symmetric  coefficients  =  r  "7,  (  are  prescribed  on  B.  Thus,  unless  such  symmetry  conditions 

hold,  the  analog  of  Piola’s  identity  (97)  does  not  necessarily  apply  for  the  divergence  of  JF~  M. 

Let  the  vector  field  A  =  AaGa  be  the  Piola  transform  of  a  =  daga: 

Aa  =JF~lAMaa.  (244) 


Taking  the  divergence  of  (244)  and  applying  the  product  rule  for  covariant  differentiation  with  definition  (217) 
results  in 

AfA  =  Aaa  =  (J~F-}AJ,A~aa  +  JF~lAMa“A  =  J[J-\jF-XAJ:Aaa  +  5“J  =  J  Vff5“,  (245) 

where  the  anholonomic  covariant  derivative  is  defined  as  [14,  19] 

VaQ  =  {-)-.aF-xam  +  {■)J-\JF-uj:A  =  (■):«  +  (-)J-l^xeABCF^B:c.  (246) 

The  second  term  on  the  right-hand  side  of  (246)  vanishes  when  the  right-hand  side  of  (243)  vanishes,  in  which 
case  (245)  becomes  similar  to  (99). 

As  a  second  example,  consider  the  following  Piola  transformation  between  contravariant  second-order  tensor 
A  =  AabGa  ®  G b  and  a  =  aAaGA  <g>  g„: 

Aab  =JF~wMdAa.  (247) 

Divergence  on  the  second  leg  of  A  is  computed  as 


AA-b  =  Aabb  =  (jF-]BJ,R~aAa  +JF-XBM~aA“B  =  J  V„5^ 

G  ... 4 

apu  r  p  BC 


=  J(daaAa  +  T ^  +  «  fcF-lBa~aCa)  +  eapxeDBC~aAaF^D(dcF^B  +  T"xFfc^). 


(248) 


Physically,  the  quantity  J  is  associated  with  volume  changes  resulting  from  mechanisms  associated  with 
mapping  F.  When  F  represents  plastic  deformation  from  dislocation  glide  [3,  43],  or  shear  associated  with 
twinning  [19,  20],  then  F  is  isochoric  and  J  =  1.  When  F  represents  thermal  deformation  [15],  and  when  such 
thermal  deformation  is  isotropic  (e.g.  for  cubic  crystals),  then  F  =  J1/31  =  J1/3g^  ga  ®  GA,  where  J  depends 
on  temperature  [25].  Mapping  F  exhibits  a  similar  spherical  form  when  it  represents  isotropic  deformation  from 
voids  [18]  or  point  defects  [17],  in  which  case  J  depends  on  the  volume  fraction  of  such  defects.  WhenF  depicts 
growth  of  biomaterials,  then  J  reflects  both  density  changes  and  mass  changes,  with  the  former  often  negligible 
in  the  context  of  tissue  mechanics  [25]. 


3.3.2.  Jacobian  determinant  of  F.  Definitions  and  arguments  of  the  previous  section  can  be  repeated  for  the  Jacobian 
determinant  of  F.  Differential  volume  elements  in  spatial  and  intermediate  configurations  are  related  by 

dv  =  J  dV,  (249) 

where  the  Jacobian  determinant  J[F(X,  t),g{x),g\  is 

j  =  l~eabceaPx  K'PVn,  =  JgTg  det F  =  ^det (gab)/ det (gaP)  det (FaJ.  (250) 

From  similar  arguments,  inverse  Jacobian  determinant  .7_l  [F  1  (.v,  t),g,g(x)]  =  1  /J  is 

J~X  =  ^eapxeabc  F-XaJ'~xpmbF~XxmC  =  yg/g  detF-1  =  ^det [gaP)f  det (gab)  det  (F_1“  ).  (251) 

Note  that  products  of  Jacobian  determinants  and  their  inverses  obey 

J  =  JJ,  (252) 


Downloaded  from  mms.sagepub.com  at  US  ARMY  RESEARCH  LABORATORY  on  September  17,  2012 


728 


Mathematics  and  Mechanics  of  Solids  I  7(7) 


From  identity  (51),  it  follows  that,  analogously  to  (242), 


dJ 

dK 


=  JF 


—  lot 


3  J 


-l 


■«>  dp-\  a 


Taking  the  divergence  of  the  first  of  (253)  results  in 

Q~XK)-.a  =  \(*PFabC 


+  F  :b:aF 


—  Co  eabcp-tP  p~lx 

—  eapxe  r  Mr 


(253) 


(254) 


Recall  from  (222)  that  spatial  gradient  y(F-1)  is  generally  symmetric  in  covariant  indices  only  when  both 

(1 12)  holds  and  symmetric  connection  coefficients  T^“  =  F "f,  (  are  prescribed  on  B.  Only  in  such  cases  does 

the  divergence  of  J~xF"a  generally  vanish,  such  that  an  identity  akin  to  (96)  holds. 

Let  the  vector  field  A  =  Aaga  be  the  Piola  transform  of  a  =  a01  ga: 


Aa  =J~xFa~aa. 


(255) 


Taking  the  divergence  of  (255)  and  applying  the  product  rule  for  covariant  differentiation  with  (218), 

A- a  =  A'a  =  0~lFaJ:aaa  +  =  J~'  V„5“,  (256) 

where  the  anholonomic  covariant  derivative  is  defined  analogously  to  (246): 

V„(0  =  (•)  :aK  +  {-VO-'F^Y.a  =  (■):«  +  (-V  ^abcF-XtiMF-\c.  (257) 

The  second  term  on  the  right-hand  side  of  (257)  vanishes  when  the  right-hand  side  of  (254)  vanishes,  in  which 
case  (256)  becomes  similar  to  (99). 

As  a  second  example,  consider  the  following  Piola  transformation  between  contravariant  second-order  tensor 
A  =  Aabga  ®  gh  and  a  =  a™ ga  ®  ga: 


Aab  =  j~lFbaaaci . 

Divergence  on  the  second  leg  of  A  is  computed  similarly  to  (248): 

Aab.b  =  Aabb  =  Q-lFbM):b~aaa  +J-XFbM~aaa.b  =  J~x  Va~aaa 

=  J~l(daaaci  +  r^af#  +  8r  b“Fba¥a)  +  eapxedbcaaaF~lfiJdcF-^b  +  T^1^1^). 

From  (79)  and  (96),  it  follows  that  covariant  derivative  operations  in  (246)  and  (257)  are  equivalent: 

v„(o  =  (-y.aK  +  yvQ-'Ky.a  =  +  {w~xJFyF-XAMyM  =  v„q. 


(258) 


(259) 


(260) 


When  F  physically  represents  elastic  deformation,  then  J  quantifies  elastic  volume  changes  typically  related 
to  pressure  changes.  The  physical  importance  of  identities  associated  with  Piola  transformations  such  as  those 
above  will  become  clear  in  the  context  of  the  balance  of  linear  momentum  discussed  in  Section  4.  Derivations 
of  Piola  transformations  and  related  identities  in  anholonomic  configurations  have  been  given  elsewhere  in 
coordinate-free  notation  [6,  40]. 


3.4.  Dislocation  theory 

Developments  of  Sections  3. 1-3.3  have  known  and  immediate  relevance  in  the  context  of  continuum  dislocation 
theory.  Connection  coefficients  referred  to  referential  and  spatial  coordinate  systems  are  introduced  as 


_  p~lA  3Dfra  —  _ pra  pP  r)aF~lA 
1  BC  —  r  .aaBr  .c  —  r  .cr  .Bapr  .a 


=  FybF~l:c  =  —F~  lamCF~ 1  dpF“c 


(261) 


Torsion  tensors  of  these  coefficients  are  related  to  the  equivalent  anholonomic  object  of  (176)  and  (184)  as 
follows: 


rj~’**A  _  T"' ../(  _  rp— 

1  BC  ~  1  [ BC\  ~  r 


Fp  Fx 

Mr  .Br  .CK  fix’ 


T£  =  rm  =  FaJr 


\.F~ 


■cKPr 


(262) 
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Riemann-Christoffel  curvature  tensors  formed  from  connections  in  (261),  which  are  said  to  be  integrable,  van¬ 
ish  identically  at  all  points  X  e  B0  or  x  e  B  (see  [27]).  Let  C  be  a  closed  curve  enclosing  simply  connected  area 
S  in  configuration  B0,  and  let  c  be  the  image  of  this  curve  in  B.  Assigning  constant  basis  vectors  ga  over  C  and 
S,  a  measure  of  incompatibility  associated  with  F  in  domain  S  is  obtained  using  (226)  and  Stokes’s  theorem: 

( tdxa  =  (f  F“a  dXA  =  f  eABCF“4.cNBdS  =  f  eABC dAF“BNc dS  =  f  aaCNcdS,  (263) 

J  c  J  C  J S  J  S  J  S 

where  Nc{X)  are  components  of  the  unit  normal  to  S.  Analogously,  letting  c  be  a  closed  curve  enclosing  simply 
connected  area  s  in  the  current  configuration, 


l 


^abcjp—  la 


crib  ds  = 


-l 


^abc  q  T7—  Iff 

€  oar 


bnc 


ds  =  Jc 


• ds , 


(264) 


with  nc(x)  the  unit  normal  to  s.  Let  F  represent  elastic  deformation,  and  let  F  represent  plastic  deformation  from 
dislocation  motion.  Dislocation  density  tensors  ceaA  and  aaa  are  related  to  torsion  tensors  of  (262)  as 

~«C  =  eABCdApaB  =  eABCp«Df..D'  ac  =  ^bc^p- 1^  =  ^bcp-laj.J .  (265) 

It  follows  from  identities  ic'p“  =  k""  of  (176),  (252),  (262),  and  Jeahc  =  eABC F°4FjFcmC  that  dislocation  density 
tensors  are  related  by 

J-'F«AapA  =JF-'amaaPa.  (266) 

Integrals  (263)  and  (264)  can  be  interpreted  as  total  Burgers  vectors  of  all  dislocation  lines  piercing  areas  S  and 
s,  respectively  [3].  Many  other  dislocation  density  tensors  can  be  derived  by  mapping  those  in  (265)  to  various 
configurations  [6,  32,  40,  42];  sign  conventions  can  also  vary  among  definitions  of  the  total  Burgers  vector 
and  dislocation  density  tensor.  A  relationship  between  torsion  and  dislocations  has  also  been  described  in  the 
context  of  gauge  theory  [44] .  In  constitutive  theories  of  crystalline  solids  of  the  gradient  type,  geometrically 
necessary  dislocation  density  tensors  such  as  those  listed  here  may  directly  affect  stored  energy  [42,  43,  45]  and 
strain  hardening  [9,  42,  43,  45]. 


4.  Balance  of  linear  momentum 

For  a  hyperelastic-plastic  material  with  uniform  properties,  it  is  shown  that  the  static  local  balance  of  lin¬ 
ear  momentum,  when  mapped  to  the  intermediate  configuration,  can  be  written  directly  in  terms  of  F  and  its 
covariant  derivatives.  Acceleration  and  body  forces  are  excluded;  these  can  be  incorporated  into  subsequent 
developments  without  difficulty. 


4.1.  Local  momentum  balance 

Let  cr(x,  t)  =  crabga  ®  g/,  denote  the  symmetric  Cauchy  stress  referred  to  configuration  B.  In  the  absence  of  body 
forces  and  acceleration,  the  local  balance  of  linear  momentum  is,  in  possibly  curvilinear  spatial  coordinates  x°, 
[37-39] 

oab-b  =  0,  (267) 

recalling  that  the  subscripted  semicolon  denotes  covariant  differentiation  as  in  (55).  Defining  the  Piola  transform 
with  respect  F  to  similarly  to  (258), 

aab  =j-lFbaPaa,  Paa=JF-\oah,  (268) 


and  appealing  to  (259),  momentum  balance  (267)  can  be  written,  upon  multiplication  by  J,  as 
0  =  Jcrabb  =  VaPaa 
=  3  aP 


jua  .  p  ..a  pafi  .  g  ..a  pb 
'  1  apr  f  p  bcr  .a 


,pca  +  jeapxedbcPaaF~lfi'd(dcF 


-lx  +  r..xp-up-u 


.b 


:*)• 


(269) 


The  first  term  in  parentheses  on  the  right-hand  side  of  (269)  is  related  to  a  dislocation  density  in  (265)  via 
chcd  d  i,  F~ 1  =  a/d .  While  the  balance  of  linear  momentum  in  anholonomic  space  has  been  considered  elsewhere 

[6,  7,  23,  40],  several  particular  forms  derived  in  curvilinear  coordinates  in  what  follows  are  believed  to  be  new. 
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4.2.  Hyperelastic-plastic  material 

Consider  a  crystalline  solid.  Let  F  physically  represent  elastic  lattice  deformation,  and  let  F  physically  repre¬ 
sent  plastic  deformation,  e.g.  resulting  from  glide  of  crystal  dislocations.  Let  'L  denote  the  strain  energy  per 
unit  volume  of  material  in  intermediate  configuration  B ,  which  is  assumed  to  be  a  function  of  the  following 
arguments: 

4/  =  vp(C,g)  =  'F[C(F,g),g],  (270) 

where  g  denotes  the  spatial  metric  tensor  of  (11)  with  components  ga/,(x),  g  denotes  the  intermediate  metric 
tensor  of  (122)  with  components  gap ,  and  where  the  symmetric  elastic  deformation  tensor 

c  =  Capga  ®  g^  =  KgahF'yga  ®  gP,  Cap  =  Fc‘agabFbmp  =  Cpa.  (271) 

The  dependence  of  strain  energy  on  other  terms  (e.g.  internal  state  variables  representing  contributions  from  lat¬ 
tice  defects,  orX  representing  heterogeneous  material  properties)  can  be  incorporated  without  severely  affecting 
subsequent  arguments.  Furthermore,  strain  energy  per  unit  reference  volume  J  T  could  be  used  alternatively 
without  conceptual  difficulties;  in  the  usual  case  of  isochoric  plastic  deformation  from  slip,  the  distinction  is 
irrelevant  because  J  —  1.  For  a  hyperelastic  response,  the  stress  obeys 


P-a=gabPha  = 


34> 


3^  3 C, 


Px 


3 FaM  dCpx  3 F“m 


dW  - 
=  2  gah-=—F' 


b 


dC, 


afi 


where  the  chain  rule  and  the  following  identity  have  been  used: 


(272) 


dCap 


Igab* 


(273) 


Constitutive  equation  (272),  in  conjunction  with  definition  (268),  is  standard  for  crystalline  solids  and  can 
be  derived  using  thermodynamic  arguments  [14,  19,  25,  33].  Functional  forms  of  T  for  nonlinear  elastic¬ 
ity  of  anisotropic  single  crystals  belonging  to  various  crystal  classes  are  available  in  the  literature  [14,  32], 
Usually,  dependence  on  g  is  not  written  explicitly;  however,  as  will  be  demonstrated  shortly,  inverse  met¬ 
ric  components  gap  are  needed  for  construction  of  scalar  invariants  of  covariant  deformation  tensor  Cap 
entering  the  strain  energy  function,  e.g.  the  trace  trC  =  C“  =  Capg and  Jacobian  determinant  j  = 
det(FaJ[det(gafe)det(^)]1/2  =  [det(C“  )]’/2. 

For  illustrative  purposes,  consider  an  isotropic  (poly)crystalline  solid  whose  nonlinear  hyperelastic  response 
can  be  described  via  the  neo-Hookean  strain  energy  potential  [46] : 

vF  =  ^[trC  +  (In/)2  -  3]  -  n  In  J  =  ^( Cap~gaP  -  3)  +  In  J  In  j  -  /r)  ,  (274) 


where  A  and  /i  are  elastic  constants.  Using  (253),  the  stress  of  (272)  for  a  neo-Hookean  material  described  by 
(274)  is 

Paa  =  /j.(Faa  -  F~laa)  +  A  F-loia  In  J  =  ti.(F“pgaP  -  F~\gab)  +  A  F~\gab  InJ.  (275) 

The  partial  derivative  of  the  stress  entering  (269)  is  then 


3 aPaa  =  PbMdbPaa  =  BPbJbFaa  -  (B  -  A  \nj)FbadbF-laa  +  U^g^dbJ.  (276) 


The  second  and  third  terms  on  the  right-hand  side  of  (276)  can  be  written  in  terms  of  the  partial  derivative  of  the 
elastic  deformation  gradient  rather  than  its  inverse  and  determinant  upon  use  of  the  chain  rule  and  the  identities 


3  F 


' —  lQf 


dF, 


_  T7 —  I**  77 

F  ./A  .cr 


d  aJ  = 


3  J 

w 


daFb  + 


3  J 


3 aig/g)  =  JP  \daF%  +  /  3fl[ln(g/g)],  (277) 


3  (g/g) 


leading  to 

3 aPaa  =  liFbadbFaa  +  (/x  -  A  \nj)(F-Uaa  dhFha  +  F-'aaFhf>dhgap  -  dbgab)  +  Xgab {F~x“c dhFcM  +  3*[ln(g/g)]}. 


(278) 


Downloaded  from  mms.sagepub.com  at  US  ARMY  RESEARCH  LABORATORY  on  September  17,  2012 


Clayton 


731 


With  stress  components  computed  generically  using  (272)  or  specifically  (e.g.  (275)),  all  that  remains  in 
(269)  is  specification  of  intermediate  connection  coefficients  T£“,  which  can  be  defined  as  discussed  in  Section 
3.2.2.  Explicit  forms  are  considered  in  what  follows. 


4.2. 1 .  Holonomic  elastic  deformation.  When  (108)  applies,  such  that  B  can  be  regarded  as  a  Euclidean  space,  then 
(190)  can  be  used  to  define  the  intermediate  connection  coefficients.  The  right- most  term  of  (269)  vanishes  from 
the  integrability  of  F~lofa  =  daxa  and  the  symmetry  of  r""  in  (190),  leaving 

3  aPaa  +  T  $P“P  +  f  tc~Fb.apCa  =  0.  (279) 


4.2.2.  Coincident  referential  and  intermediate  coordinate  frames.  Next  consider  the  case  when  (109)  may  apply,  such 
that  single-valued  coordinates  xa(X,  t )  continuously  differentiable  with  respect  to  XA  may  not  exist  in  B.  Let 
identical  coordinate  systems  and  metric  tensors  be  assigned  at  to  each  material  particle  X  in  configurations  B0 
and  B  [36],  such  that  the  appropriate  connection  coefficients  are  given  by  (196).  In  this  case,  linear  momentum 
balance  (269)  becomes 


daPaa  +  BC~F^B.aPaP  +  f 

+  (detF  )eanedbcPaaF~X\ 


>apb  pc i 


be 


(3  CF 


'  _i_  f  X  ?CG 

+  dA  X  r  rc 


■AF-'B_CF-Umh)  =  0. 


(280) 


4.2.3.  Coincident  spatial  and  intermediate  coordinate  frames.  Now  consider  the  case  when  (113)  may  apply,  such  that 
single-valued  coordinates  xa(x,  t )  continuously  differentiable  with  respect  to  xa  may  not  exist  in  B.  Identical 
coordinate  systems  and  metric  tensors  are  assigned  to  each  spatial  point  x  in  configurations  B  and  B,  so  that 
appropriate  connection  coefficients  on  B  are  given  by  (206)  (see  also  the  Appendix  for  a  specific  example  of 
cylindrical  coordinates).  Linear  momentum  balance  (269)  becomes 

3 aPaa  +  f  bcF^S^  +  Pca) 

+  (detF )eafixedbcPaaF-lliJdcF-^b  +  $*5*  ^F~Umb)  =  0.  1 1 

Choice  (206)  appears  favorable  over  (196)  in  the  context  of  the  balance  of  linear  momentum  because  (281) 
involves  only  one  set  of  connection  coefficients  and  one  deformation  map  (F),  in  contrast  to  (280)  which 
involves  two  different  sets  of  connection  coefficients  and  both  deformation  maps  F  and  F.  When  (278)  applies, 
for  example,  momentum  balance  (281)  can  be  expressed  completely  in  terms  of  the  elastic  constants,  elas¬ 
tic  deformation  F  (including  its  gradient,  inverse,  and  determinant),  spatial  metric  g,  and  spatial  connection 
coefficients  p  the  latter  which  are  obtained  from  partial  derivatives  of  g  through  (53). 


4.2.4.  Cartesian  intermediate  coordinates.  Finally  consider  the  case  wherein  Cartesian  basis  vectors  are  used  on  B, 
so  that  (214)  holds  meaning  Tf  =  0.  The  linear  momentum  balance  (269)  reduces  to 

3 aPaa  +  f  bacFhMPca  +  (detF )e0lPxedbcPaaF-lfimddcF~lx.b  =  0.  (282) 

When  Cartesian  spatial  coordinates  are  used  and  when  F-1  is  integrable  as  in  (279),  then  (282)  reduces  to 
daPaa  =  0. 

5.  Conclusions 

The  possibly  anholonomic  intermediate  configuration  of  a  deformable  body  has  been  examined  from  the  per¬ 
spectives  of  tensor  calculus  and  differential  geometry.  Various  choices  of  extrinsic  coordinate  systems  for  the 
intermediate  configuration  have  been  considered;  corresponding  metric  tensors,  connection  coefficients,  torsion, 
and  curvature  have  been  derived.  Partial  and  total  covariant  differentiation  with  respect  to  anholonomic  coordi¬ 
nates  have  been  defined.  It  has  been  shown  that  when  the  same  curvilinear  coordinate  systems  are  prescribed 
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on  intermediate  and  Euclidean  reference  or  current  configurations,  the  torsion  and  curvature  of  the  connection 
coefficients  in  the  intermediate  configuration,  which  are  generally  time  dependent,  need  not  vanish.  Conven¬ 
tional  identities  that  hold  in  reference  and  current  configurations  such  as  vanishing  divergence  of  the  curl  of 
a  vector  field  and  vanishing  curl  of  the  gradient  of  a  scalar  field  do  not  necessarily  apply  in  the  intermediate 
configuration.  The  balance  of  linear  momentum  has  been  derived  for  a  hyperelastic-plastic  material  in  general 
curvilinear  coordinates.  For  neo-Hookean  elasticity,  the  choice  of  coincident  spatial  and  intermediate  coordinate 
systems  enables  the  stress  divergence  to  be  expressed  completely  in  terms  of  elastic  constants,  spatial  gradients 
of  elastic  deformation,  and  the  spatial  metric  tensor  and  its  spatial  derivatives. 
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Appendix:  Anholonomic  cylindrical  coordinates 

As  an  illustrative  example  of  curvilinear  intermediate  (anholonomic)  coordinates,  consider  cylindrical  spatial 
coordinates  mapped  to  the  intermediate  configuration.  Because  cylindrical  spatial  coordinates  exhibit  rela¬ 
tively  simple  forms  for  metric  tensors  and  connection  coefficients,  corresponding  quantities  mapped  to  the 
intermediate  configuration  can  be  derived  by  inspection. 

In  this  example,  (146)— (15 1)  apply  for  basis  vectors  and  metric  tensors.  Since  the  angular  coordinate  in  cur¬ 
rent  configuration  B  is  labeled  9,  overbars  are  used  here  to  denote  particular  (anholonomic)  coordinates  referred 
to  the  intermediate  configuration  B.  Specifically,  holonomic  coordinates  on  B  are  denoted  by  the  usual  (r,  9,z), 
while  anholonomic  coordinates  on  B  are  denoted  by  (r,  0,  z).  In  what  follows,  indices  (r,  9,z )  and  (r,  6,  z)  refer 
to  specific  coordinates  and  are  exempt  from  the  summation  convention.  Generic  indices  ( a ,  b,c,d ...)  and 
(a,  /J,  x,  S  . . .)  refer  to  free  spatial  and  intermediate  quantities,  respectively,  and  are  subject  to  summation  over 
repeated  indices. 


A.I.  Spatial  components 

Quantities  associated  with  the  usual  holonomic  spatial  coordinates  are  considered  first  [39],  In  three  spatial 
dimensions,  cylindrical  coordinates  are 

(x1,x2,x7>)  — »•  (r,  9,z),  r  >  0,  9  e  (— n,n].  (283) 
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The  squared  length  of  a  differential  line  element  dx  is 

dx  -  dx  =  (dr)1  +  ( rdd )2  +  ( dz )2. 

The  metric  tensor  and  its  inverse  are 
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Determinants  of  the  metric  tensor  and  its  inverse  are 

g  =  det(gafe)  =  r2,  g~l  =  det(gfl*)  =  \/r2. 


(284) 


(285) 


(286) 


A  differential  volume  element  is 

dv  =  ^/g  dx 1  dx 2  dx 3  =  rdrdO  dz.  (287) 

Connection  coefficients  from  (53)  are 

f  "0  =  f  Or  =  l/r>  f  00  =  ~V’  f  be  =  0  0therwise-  (288) 

It  is  remarked  that  connection  coefficients  (288)  are  symmetric  in  covariant  indices  (i.e.  null  torsion),  and  that 
the  Riemann-Christoffel  curvature  tensor  formed  from  these  coefficients  vanishes.  Let  x  denote  the  position 
vector  in  Euclidean  space.  Natural  basis  vectors  are 


gr  =  drx,  ge  =  3 ex,  g_-  -  3zx. 


(289) 


A.2.  Intermediate  components 

Intermediate  basis  vectors,  from  (289)  and  (146),  are 

g  r(x)  =  9,-X,  gg(x)  =  dgX,  gz(x)  =  3ZX. 
The  metric  tensor  and  its  inverse,  from  (285)  and  (147),  are 
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\/r2 
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0 

1 

(290) 


(291) 


Determinants  of  the  intermediate  metric  tensor  and  its  inverse,  from  (286)  and  (149),  are 

g  =  r2,  g_1  =  \/r2.  (292) 

To  second  order  in  dx,  from  (226),  the  squared  length  of  an  intermediate  line  element  is 

dx  •  dx  —  (F~laml. dr  +  F~la£  dd  +  E’“1“  dz)gap(F~lpml. dr  +  F~l d6  +  F~lfl dz).  (293) 

An  intermediate  volume  element  is,  from  (151),  (249),  (251),  and  (287), 

dV  =  (detf~l)rdrdO  dz.  (294) 

In  the  present  example,  definitions  (206),  (208),  and  (209)  apply  for  connection  coefficients,  torsion,  and  cur¬ 
vature  referred  to  the  intermediate  configuration.  From  (288)  and  (206),  anholonomic  connection  coefficients 
consist  of  up  to  nine  non-zero  components: 


F-/r  T--=F-/r  T“-  =  F- /r 

1  rt'  ’  1  QQ  1  Ql  '  ’  1  1  z!  '  ■> 


zO 


(295) 
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T-  =  F-  /r  T-  =  F  -  /r  T"  =  F-  /r 

T-f  =  —F-r  f-  =  -F°r  F"'  =  -F°r 

w  ’  oe  1  o  '>  ~e  z 

Possibly  non-zero  skew  covariant  components  of  these  coefficients  are 

^  ( T?r  T?8\  t -’■.0  ^  T7r . 

rpe]  _  2r  F§’’  ri^l  “  2r’ 

F"0  — 

1  k'']  2r  z  ’ 

p..r  _ _ —F®  p-p  _ _ —F® 

m  -  2  PI-  2 

Torsion  can  then  be  computed  from  (208),  noting  that  the  anholonomic  object 

Spatial  connection  coefficients  (288)  have  the  following  non-zero  skew  partial  derivatives: 


a  g  ..8  _ _ J_  a  g  ..r 

a[rrey—  2  r2’  ['T0F 


1 

2 


(296) 

(297) 

(298) 

(299) 

(300) 

(301) 

(302) 

(303) 


Therefore,  possibly  non- vanishing  components  of  the  anholonomic  curvature  (209)  follow  as 

2 

R—e  —  _ pr  pe  —  iFr  F6 

KPxr  ~  r2r-[ F-xY  Kpxe  -  zr-[F-xY 

Written  out  completely  with  no  free  indices,  (303)  yields  the  components 


n-.T-  _  pr  pfi  _  pr_p0_  n...r  _  pr_p8_  _  pr  p0  n...r  _  pr_p8_  _  pr  p8 

?88  -r  .0  .9  -r’  8~z8  .8  -z  1  -z‘  . 8 5  zr8  ■-  -r  -r  .z' 


(305) 


Clearly,  skew  components  of  the  connection  (and,  hence,  the  torsion)  as  well  as  the  curvature  need  not  vanish 
even  if  the  anholonomic  object  of  (301)  is  identically  zero. 
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